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V. PHYSICAL AND MATHEMATICAL MODELS 


5.1 SINGLE-PHASE FLOW MODELS 

In this section, the equations which govern the distribution of the mean quantities are summarized. 
These equations are derived from the conservation laws of mass and momentum using time averaging 
and are expressed in tensor notation for steady and constant density flow as 



dXj 


ftUjUj) i ap 

ax, p ax* 


3X; 


(u iUj ) 


( 6 ) 


where Uj and uj are the mean and fluctuating velocities along the Xj direction, respectively, p is the mean 
pressure, and the bar is used to denote time-averaged quantities. 


As a consequence of the nonlinearity of Equation 6, the averaging process used introduces unknown cor- 
relations UjU, which can be made known through the assumption of turbulence modeling. 


Three different types of turbulence closures are investigated, namely, the standard k-e model, algebraic 
second-moment closure (ASM), and differential second-moment closure (DSM). 

The k-e model is a simple closure based on the gradient transport relations. In this model, the turbulent 
fluxes are related to the mean fields through the assumption of an isotropic eddy viscosity and a turbu- 
lent Prandtl/Schmidt number as 


-pUiU,=p, 


3Uj dU, 



(7) 


The eddy viscosity (p^) is obtained from the turbulent kinetic energy (k) and its dissipation rate (e) using 
the relation 


p t =c^pk 2 /e (8) 

In order to close the set of Equations 5 through 7, two additional equations governing the transport of k 
and e are required. These are 
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where ojc and o e are turbulent Schmidt numbers and C E l and Q2 are model constants. The constants 
used in this model have been taken from Launder and Spalding (1974*) and are given in Table 5.1-1. 


Table 5.1-1. 

Values of constants in the k-e model. 

-Cp- -Cel- -Ce2- -Qk_ 


0.09 1.44 1.92 1.0 1.3 


The k-e model is the simplest model which is suitable for recirculating flow calculations. It allows the 
characteristic length scale of a wide range of complex flow fields to be determined. The k-e model has 
been used with success in the calculation of various free shear flows and recirculating flows with and 
without swirl (e.g., Rodi, 1980). However, in flows with significant streamline curvature, the isotropic 
eddy viscosity assumption may not be able to describe the turbulent diffusion effects adequately. The 
axisymmetric form of the turbulent flow equations is given in Appendix A for the k-e model. 

To allow for the nonisotropic behavior of the eddy viscosity and to account for the effect of body forces 
(e.g., buoyancy, rotation), the k-e model is refined by introducing ASM. This model is based on a simpli- 
fication of the Reynolds stress transport equation which relates the individual stresses to mean velocity 
gradient, turbulent kinetic energy, and its dissipation rate by way of algebraic expressions. The ASM 
model adopted here is based on Rodi’s hypothesis (Rodi, 1976) which approximates the convection and 
diffusion transport of turbulent stresses in terms of the transport of k. 

The result can be summarized as 



( 2 \ 

1-C 2 


7! 7~K 

pe 
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( 11 ) 


where 


dUj dU: 

,i=pUiUk axi'" pUiUk ax7 (12) 

p k=^Pii (13) 

Cl and C2 are model constants and ajj is the nondimensional measure of anisotropy and is given by the 
following expression 


a ii = 



(14) 


* References for Section V are listed at the end of the section. 
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Since the quantities k and e arc present in these equations, their transport equations also have to be 
solved. These are determined by 


IT dk d 
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’ dX dX 


C w p— u i u n 


dk 




i. d£ d 

pUj = 

1 3Xj axj 


CcP-«iUr 


ax, 

ae 


-P k -pe 


' c e 1 "ax 


n ) 


+ C tl ^Pk-C, 2 p^— 
k k 


Ck, C zj Qi, and Ce2 are all model constants and are given in Table 5.1>1I. 


(15) 


(16) 


Table 5.1 -II. 

Values of constants in the ASM and PSM closures. 


Cel 

1.44 

C E 2 

1.92 

C E 

0.18 

Ck 

0.22 

Cl 

1.8 

C 2 

0.6 

Clw 

0.05 

C2w 

0.006 


The k-e and ASM models assume that the local state of turbulence can be characterized by one velocity 
scale. In order to allow for the different development of the various Reynolds stresses representing vari- 
ous velocity scales in complex flows and to account properly for their transport, models which employ 
transport equations for the individual stresses must be applied. 

The Reynolds stress equations can be written in tensor notation form as 

P U k “ dij = P«j + ^ij “ P € ij < 17 > 

dX k 

Here, Pjj is the production of Reynolds stress ujuj, ejj represents viscous dissipation, <(>ij controls the redis- 
tribution of turbulence energy among the normal stresses through the interaction of pressure and strain, 
and djj stands for turbulence diffusion. Since Pjj is exact, it does not need modeling. However, closure 
assumptions are required for dij, 4>ij, and ejj. The assumption of local isotropy allows the dissipation ten- 
sor to be approximated by 

t,=|s i|£ (18) 

where £ is the turbulence energy dissipation rate. The diffusion term is approximated by the gradient- 
diffusion model of Daly and Harlow (1970) 
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(19) 
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where Ck is a model constant, and k = 1 /2 u.u, is the turbulent kinetic energy. The pressure redistribu- 
tion term (<|>ij) is modeled in three parts: fojl resulting from purely turbulence interactions known as "re- 
tum-to-isotropy," 4>ij2 involving interactions between the mean strain rate and turbulence known as 
"rapid" part, and foj w representing the effects of rigid boundaries on both foji and fop. The presence of a 
rigid wall affects the flow field near that region by impeding the transfer of turbulent energy from the 
streamwise direction to that normal to the wall, and as a result reduces the relative magnitude of the 
shear stress. In the present study, Rotta's linear model (Rotta, 1951) for the turbulence part of foj is 
adopted ’ 


♦iji = ~ c iP- 


u i u r -5 j j k 


(20) 


The simple linear form for foji is widely accepted and used dispite the fact that the actual "return-to- 
isotropy" process is highly nonlinear (Bradshaw, 1968). More sophisticated nonlinear forms, such as 
Lumely and Khayeh-Nouri's proposal (1974) have been suggested, but these have shown no significant 
improvement over Rotta's proposal. The rapid part is approximated using the simple model suggested 
by Naot et al. (1970), known as the isotropization production (IP) model ’ 


♦hi— C, P u —8,^ 
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( 21 ) 


A more sophisticated version of foj2 is the linear quasi-isotropic (QI) model (Launder et al, 1975) 
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This model includes both the symmetric and antisymmetric mean strain effects on redistribution model- 
ing. 


The effects of solid boundaries on pressure redistribution term are included using the wall correction 
proposed by Launder et al. (1975), or 
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kl 


U i U r 3 5 M P k 



(24) 


where X n is the normal distance from the wall and the model constants Cl w and C2w are specified in 
Table 5.1 -II. The modeled Reynolds stress transport equations in axisymmctric coordinates (x, r) are 
given in Appendix B. 
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5.2 TWO-PHASE FLOW MODELS 

The mathematical formulations for the particle/ flow problems consist of Eulerian conservation equations 
to the continuous gas phase and Lagrangian equations to discrete particle motion. The continuous gas 
phase is coupled to the discrete phase by incorporating additional source term. 

5.2.1 Dispersed Phase 

The model employed for the dispersed phase calculations is based on the stoichastic Lagrangian formula- 
tion reported by Gosman and Ioannides (1983). In this approach, effects of the gas phase turbulent fluc- 
tuations on particle drag and dispersion are considered. Because of the large number of actual particles 
in the spray, the dispersed phase is characterized by individual spherical particles, termed "computa- 
tional particles." Each of these computational particles represents a group of particles all having the same 
initial size, velocity, and temperature. The commonly adopted equation of motion of n-th computational 
particle in dilute spray, assuming large particle-to-fluid density ratio, is given by 

dv L= (u iZ vr) +g „ 

dt xd 

where Uj and Vj are the ith components of instantaneous gas and particle velocity, respectively, g is the 
gravitational acceleration, and td * s dynamic relaxation time of particle defined as 

T .= — 4I ?P P P (26) 

3C D p|U-V n | 

Here, Dp and pp are particle diameter and density, respectively, Cd is the drag coefficient, p is the fluid 
density 


C D = (24/Re)(l + 0.1315[Re] 0 82 ' 0 05w ), 0.01 < Re 5 20 (27) 

Cq = (24/Re)(l + 0.1935lRe] a6305 ), 20 < Re < 260 (28) 

where w = logioRe and the particle Reynolds number (Re) is defined as 


p|u-V n D 

Re = J 

P 


(29) 


where p is the fluid density. Integration of acceleration from Equation 25 results in the velocity compo- 
nents of the particle. The position of each particle group can be found by integrating the equation 


(30) 

dt 1 

where Xj is the particle position vector. A particle is assumed to interact with an eddy for a characteristic 
time x. The interaction time is determined by the minimum of either the residence time (tj) which is the 
time required for a particle to cross an eddy or the eddy lifetime (te) which is the time that particle re- 
mains within the eddy during the whole of its lifetime. The time scales are estimated with the assump- 
tion that the characteristic size of an eddy is the dissipation length scale L^, given by 
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(31) 


Ue-c^V^/e 
te and t r are then estimated as 

t e =Lg /V(2k/3) (32) 

t r = L e /|u-V n | (33) 

hence, - - - 

t - irdnde, t r ) “ ' (34) 


The velocity fluctuations associated with each eddy are found by making a random selection from the 
probability density function of velocity. The turbulence is assumed to be anisotropic if the DSM model is 
applied with fluctuating components having a Gaussian distribution. The standard deviation for each 
distribution is taken to be its respective root mean square (rms) velocity fluctuation obtained from the 
Reynolds equations. The distribution is randomly sampled when a particle enters the gas velocity field to • 
obtain the instantaneous gas velocity. However, in the k-e model, the rms velocity fluctuation is chosen 
randomly from an assumed isotropic Gaussian distribution with mean square deviation 2/3k. 

5.2.2 Continuous Phase 

For high Reynolds number turbulent flow, neglecting the volume fraction of the particle phase, the time 
averaged mass and momentum equations can be expressed in Cartesian tensor notation as 

4-<P u i> = ° (35) 

dXj 



( p U,U i )-^--^-(pu,u j )-S, 


ax, 


(36) 


where Spi is the interfacial drag force resulting from interaction between particles and the carrier phase. 
The interfadal drag force per unit volume is given by 


Sp< = 2^ ct n c D “~|u - v n (u, - vr ) 


(37) 


where a n is the volume fraction of particle group n that passes through the computational cell and is 
given by 



Here/ N is the number of particles represented by the trajectory n, and AV is the computational cell vol- 
ume. 
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Two different types of turbulence closures are investigated. One is a simple closure based on the gradient 
transport relation (Equation 7). Two additional partial differential equations are also needed to obtain k 
and e (Mostafa and Mongia, 1988). These are 
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F n and tl are the interphase friction coefficient and the continuous phase Lagrangian time scale, respec- 
tively, given by 



x L = 0.35^ (44) 

L t 


Another type of closure is based on solving additional transport equations for UjUj without invoking the 

gradient transport assumption. The Reynolds stress equation suggested for two-phase flow can be writ- 
ten in tensor notation as 



= P ii+4»i j -pe i j + P i j 


(45) 


where Pjj, eij, and djj represent production, dissipation, and diffusion of Reynolds stresses (Equations 12, 
18, and 19). Pjj is an extra term that accounts for the generation of turbulence due to the particles motion. 
The model for Pjj tensor is suggested as (Lopez et al, 1990) 


- 4/500 

p ij = 0 3/5 0 
L 0 0 3/5J 


(46) 


where Pj is the turbulence production due to the particles work as they move through the continuous 
phase and is given by 



(47) 
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where C3 = 0.02 is a model constant indicating only 2% of the turbulence created by particles is trans- 
ferred into the large eddies. 

The pressure redistribution, is modeled in three parts. These are: <&iji resulting from purely turbu- 
lence interaction, <&ij2 involving interactions between the mean strain rate and turbulence, and 4>\n, as 
suggested by Lopez et al (1990), representing the effects of particles on both and f 

♦ ir -C,pt(5^.2 {lj |,).c 2 (P 1| -|6 ij P k )-C 2 (P r 26 ij P k ) (48) 

The additional unknown, e, is determined by the following equation 



(49) 
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53 MATHEMATICAL FORMULATION 


The calculation procedure used in this study is based on the primitive variable formulation of the Navier- 
Stokes equations. The conservation equations are discretized using a control-volume approach. The 
coupling between the continuity and momentum equations is handled via the SIMPLER algorithm. The 
procedure is described in detail in Patankar (1980), Varejao (1979), and Karki et al (1988). 

The conservation equations for all dependent variables may be expressed in the following general form 


A (pW .A, p v., 


_a 
ax i ax 


r* 




+ s 


(50) 


where $ is the particular variable of interest, r is the diffusion coefficient, and S is the source term. 


53.1 Discretization 


Equation 50 can be written as 


a t a , 

w Jx+ aY , y i 


(51) 


where J x and Jy are the total (convection and diffusion) fluxes defined by 



(52a) 

Jy'PVf-r^ 

(52b) 


The integration of Equation 51 over the control volume surrounding the grid point P (Figure 5.3.1-1*) 
gives 

<Jx,e * Jx,w> AY + (J y , n - Jy ;S > AX = SAX AY (53) 

A discretization scheme is needed to relate the flux at each control-volume face to the values of the de- 
pendent variable at the neighboring grid points. The results presented in this report have been obtained 
using the power-law differencing scheme and flux-spline scheme. A brief description of these schemes is 
presented next. 

53.2 Power-Law Differencing Scheme 

This scheme is based on a curve fit to the exact solution of the one-dimensional convection-diffusion 
equation without a source. Since this formulation is based on a purely one-dimensional flux balance, it 
leads to significant numerical errors in the presence either of strong source terms, or of crossflow gradi- 
ents in multidimensional flows coupled with the grid-to-flow skewness. The flux-spline scheme includes 
these effects in the interpolation profile between the grid points. 


* Figures for Section V appear at the end of each subsection. The figure number identifies the subsection 
in which the figure is discussed. 
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The flux-spline scheme considered here is based on the assumption that within a control volume the total 
flux in a given direction varies linearly along the coordinate direction. For example, the flux in the x-di- 
rection for the control volume around the grid point P (Figure 5.3.3-1) is given by s “ 

The integration of Equation 54 leads to the following expression for the variation of ♦ in the x-direction 

♦ *a + bexpf)UX/D + c> ($ 5 ) 

where the constants a, b, and c for a given control volume can be expressed in terms of J x>e , J X/W , and $p. 

Equation 55 gives the variation of $ within a control volume. For two adjacent control volumes the $- 
profiles are such that they imply the same total flux at the common interface. In addition, these profiles 
must also give a unique value of $ at the common interface. This continuity-of-^ (spline continuity) con- 
dition for the interface between the grid points W and P can be expressed as 

Jx,w = ((^x,w$ W * ^w^p) + ®x,w (Jx,w * Jx,e^ + (--x,w (J X/ w * Jx,ww) (56) 

Here, the expression (D X/ w4>W " Ex,w4>p) is identical to that obtained from the lower-order exponential 
scheme (e.g., Patankar, 1980) which is based on the assumption that the total flux is uniform within a con- 
trol volume. The extra terms involving Bx and C x result from the linear variation of flux. For ease of pre- 
sentation, Equation 56 is rewritten as 

A 

Jx,w = (^x,w$W ' ^x,w^P) + Jx,w (57) 

It should be noted that additional terms such as J X/W involve the difference in flux values at adjacent faces 
of the control volume. That there is a difference in flux indicates the presence of a source term and/or 
multidimensionality (a change of flux in one direction is felt as a source term in another direction). 

Similar expressions can also be derived for fluxes in other coordinate directions. Substituting these ex- 
pressions in Equation 53 and utilizing the discrete form of the continuity equation, the following dis- 
cretization equation for <|> is obtained 

a p^p = E a nb^nb + ^ + 5 (58) 

The values of the influence coefficients a n b are identical to the coefficients obtained from the exponential 
scheme. The contribution of the flux-spline formulation is contained in the term S, which is given 

5 = Ox,w " Ix,e) A y + (Jy,s ' Jy,n) Ax (59) 

A two-dimensional situation is governed by three field variables: 4>, J X , and Jy. The three sets of equa- 
tions that determine these variables are 

(1) the conservation equation for $ 

(2) the spline-continuity condition in the x-direction 

(3) the spline-continuity condition in the y-direction 
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A. A 

The solution of these equations is obtained in an iterative manner. In the beginning, J x and J y are set 
equal to zero, then the conservation equation for $ reduces to the lower-order formulation and can be eas- 
ily solved. The solution leads to new estimates for the fluxes J x and Jy from which new J x and J y can be 
calculated. The ^-equation is now solved with the flux-spline contribution to the source term. This pro- 
cess is repeated until convergence is achieved. 
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Figure 5. 3. 3-1. One-dimensional condition. 
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VI. MODEL VALIDATION 


This section presents the results of the comparisons of the various models for single-phase and two-phase 
turbulent recirculating flows. All the experimental data have been fully presented and discussed in 
Chapter IV. Predictions obtained with each of the turbulence models are discussed and compared with 
experimental data. The goals of this study require a careful selection of the test cases. They have to pro- 
vide reliable mean flow and turbulence data in the recirculation and recovery regions. A detailed specifi- 
cation of the flow parameters in the upstream region is also essential since these are used as inlet condi- 
tions to start the computations. The test problems used range from simple flows to more complex flows 
to encompass the range of complexities involved in combustor flows. These are 

• single round jet 

• single annular jet 

• single swirling annular jet 

• coaxial jets 

• coaxial jets with swirling annular flow 

• airblast injector 

All data sets meet the criteria for acceptable measurements in terms of adequate experimental facility, ap- 
propriate instrumentation, and agreement with generally accepted flow trends. 

6.1 SINGLE ROUND JET 

This section presents the comparison of computational results with experimental data. Both the confined 
and unconfined configurations with and without the presence of glass beads are investigated. Bench- 
mark quality data were obtained by using a two-component phase/Doppler technique. The theoretical 
approach is based on a stochastic Lagrangian treatment for the discrete phase combined with an Eulerian 
description of the fluid field and using the k-e or differential second-moment (DSM) closure. 

6.1.1 Unconfined Single Round let 

The case of round jet flow configuration was selected for model validation. The glass bead injector (D = 
25.3 mm i.d.) with no inlet swirl was directed vertically downward within a 457 mm^ wire mesh screen. 
The sketch of the test section is shown in Figure 6.1. 1-1*. Data were obtained at seven axial stations 15, 25, 
35, 50, 75, 150, and 300 mm from the exit plane of the injector. In addition, measurements were made 
close to the exit plane of the rejector to provide initial conditions (Table 6.1-1). 

The carrier phase governing equations are solved numerically using the marching finite-difference solu- 
tion procedure (Spalding, 1978* ). The present calculations were obtained using a fine grid with 100 
cross-stream grid points and marching step sizes limited by 3% of the current radial grid width or an en- 
trainment increase of 3%, whichever is smaller. The ordinary differential equations governing particle 
motion are solved using a second-order finite difference algorithm. Ten thousand particles are used for 
the stochastic (ST) treatment, whereas 200 particles are computed when the deterministic (DT) method is 
compared with the ST. 

The calculations start from the first experimental location (1 mm from the nozzle exit) where measured 
mean and root mean square (rms) velocity profiles for both gas and particles are available. The inlet pro- 
file for turbulent kinetic energy dissipation rate is obtained from the measured turbulent kinetic energy, 
turbulent shear stress, and axial velocity gradient. 


* Figures for Section VI appear at the end of each subsection. The figure number identifies the subsection 
in which the figure is discussed. 

* References for Section VI are listed at the end of the section. 
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Table 6.1-1. 

Experimental flow conditions at Q.Q4P downstream of pipe exit* 


Parameter 

Gas-phase (air) 

Centerline velocity, U Zf o, m/s 
Density, n, kg/nr* 

Mass flow rate mi, kg/s 
Reynolds number Re = 4mi/pmjD 
Particle-laden (glass beads) 

Centerline velocity, V Z; o, m/s 
Density, 12 , kg/nv* 

Mass flow rate m 2 , kg/s 

Ratio to particle-to-gas mass flow rate, LR 


Single- 
phase jet 

Particle-laden jet 
Case 1 Case2 

4.74 

4.70 

4.20 

1.178 

1.178 

1.178 

0.0021 

0.0021 

0.0021 

5712 

5712 

5712 


4.41 

339 


2500.0 

2500.0 


0.00042 

0.0021 


0.2 

1.0 


* D = 0.0253. 


Predictions of a single-phase fluid flow field are first compared with data. The performance of the turbu- 
lence model in the jet developing region and the effect of the inlet dissipation rate on predicted airflow 
field are assessed. Both the standard k-e and its modified version (Spalding, 1978) for round jet are eire 
ployed. In the modified k-e model, the effect of the centerline velocity gradient on and C p? is given by 

cp = 0.09 - 0.04fj 

ce 2 = 1.92 - 0.0067fi (61) 


(60) 


where 


f, =|0.5R- 


dU c 

dU c 

[0.2 

d z 

dz 


Uc- 




(62) 


Uc and U« are the axial velocities of the fluid at the jet centerline and the ambient stream, respectively, 
and R is the local jet width. This modification allows the atainment of the correct spreading rate of a self- 
similar round jet. 

Figures 6.1. 1-2 and 6.1. 1-3 relate to the measurements of a single-phase jet and present comparisons with 
calculations using both k-e models. The results are plotted in a dimensionless form versus r/r Q , where r G 
is the radius of the nozzle. This way the jet spreading can be seen from the mean axial velocity profiles. 
The mean axial gas velocities are normalized by the centerline velocity at the inlet of the nozzle, U Z/0r so 
that the jet centerline velocity decay can be illustrated in the same figure. All other quantities are normal- 
ized by the local mean centerline velocity. Figure 6.1. 1-2 presents mean axial velocity, turbulent kinetic 
energy (k), and shear stress, while Figure 6.1. 1-3 shows the three normal stresses and k. The gas rms ve- 
locities were calculated using isotropic assumption. 

Both models yield good overall agreement with measurements dose to the nozzle exit, but for the region 
downstream, the standard model performs better. In the first region, the decay of the centerline velodty 
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is very small, therefore, the two sets of model constants are quite the same. When the axial velocity 
started to decay substantially (at about Z/D = 5) Cjj, and c^2 of the modified model decrease to values that 
give the observed spreading rate for self-similar round jets. For instance, at Z/D = 7, C|x reduces to a 
value of 0.064, whereas c^2 becomes 1.87. These low values decreased the turbulent diffusion and, hence, 
the decay of the centerline velocity. The modified model also underpredicted the level of turbulence en- 
ergy. The reason may be attributed to the inlet dissipation rate. Therefore, calculations were made with 
the standard model, but with two inlet dissipation levels, as given by Eq and 0.5eo where Co is the dissi- 
pation reference conditions obtained from the measured profiles of mean velocity, turbulent kinetic en- 
ergy, and shear stress at 1 mm downstream of the nozzle exit. As seen from Figure 6.1. 1-4, Italving the 
inlet dissipation rate improves the predictions of the kinetic energy of turbulence. 

The comparison of two versions of the k-e model indicates that the results of the standard model agree 
better with experimental data, especially for the mean axial velocity. Therefore, this model was extended 
for two-phase flow calculations. 

For particle-laden jets, two mass loading ratios (LR), defined as the ratio of particle-to-gas mass flow rate 
at the inlet plane, LR = 0.2 and 1 .0, were considered. To distinguish between the effects of mean and fluc- 
tuating gas velocity in particle transport, predictions using ST and DT treatments are compared with the 
measured data. The main difference between the two treatments is that the first considers the effect of 
gas turbulence on particle motion, whereas the second completely ignores it. 

Figure 6.1. 1-5 shows the measured and predicted (ST treatment) distributions of mean centerline veloci- 
ties and particle number density normalized by their corresponding values at the nozzle exit. Because of 
the slow decay of the particle mean velocity, there is a momentum transfer from the dispersed phase to 
the gas which causes an increase in the latter velocity compared with the corresponding single-phase val- 
ues. This change in the gas-flow properties could also be attributed to the gas turbulence attenuation 
caused by particles. Gas turbulence attenuation causes a reduction in the jet spreading rate that results in 
less decay of the gas centerline velocity. Figure 6.1 .1 -5 (a) shows a progressive increase in the gas center- 
line velocity with the particle loading ratio LR, which might be explained by the fact that both the mo- 
mentum transfer and turbulence modulation are proportional to LR. It can be seen from Figure 6.1 .1-5 (b) 
that the particle number density is decaying much faster than the particle velocity. This means that the 
spreading of the dispersed phase dilutes the particle concentration but does not necessarily decelerate the 
particles. 


Figures 6.1. 1-6, 6.1 .1-7, and 6.1. 1-8 correspond to the measurements of a particle-laden jet with LR = 0.2 
and indicate the extent to which the models described allow realistic calculations. Figure 6.1. 1-6 presents 
the mean velocities of both phases, whereas Figures 6.1 .1-7 and 6.1. 1-8 show the rms values for the parti- 
cles and gas and the shear stress of the latter. In these calculations, the turbulence model presented for 
two-phase flows was used. It can be seen from Figure 6.1 .1-6 that both ST and DT treatments yield nearly 
the same results for gas quantities. However, for particle quantities the ST provides good predictions 
compared with the experimental data, whereas the DT performs quite poorly for the particle flow proper- 
ties. According to the latter, a particle moves radially due to its initial mean radial velocity and/or the 
mean radial gas velocity, both of which are very small compared with the axial component This might 
explain the narrow distribution of particle mean axial velocity and number density predicted by the DT. 

Figure 6.1. 1-7 shows comparisons between predicted and measured values of gas turbulence kinetic en- 
ergy and shear stress. It also presents the rms axial gas velocity and its corresponding value for the parti- 
cles. If Figure 6.1. 1-7 is compared with Figure 6.1. 1-2, some reduction in the gas kinetic energy of turbu- 
lence caused by the particles could be observed. This phenomenon is more pronounced at LR = 1.0 and 
will be discussed below in connection with the results of that case. In Figure 6.1. 1-8 the comparisons are 
made for the radial and azimuthal rms velocities of both phases. It is seen that the ST predictions of rms 
particle velocities are in very good agreement with the data; v z ^ is somewhat underpredicted near the jet 
centerline and overpredicted at the outer boundary. This is in agreement with the calculations of Shuen 
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et al (1983) and Bulzan et al (1987) who studied particle-laden jets under somewhat different flow condi- 
tions. 

Figures 6.1. 1-9, 6.1.1-10, and 6.1.1-11 relate to the measurements of particle-laden jets at LR = 1.0 and pre- 
sent comparison with calculations. Figure 6.1. 1-9 shows mean axial velocity values for both phases and 
particle number density, and Figures 6.1.1-10 and 6.1.1-11 show the Reynolds stresses. It can be seen from 
these figures that the DT is inferior to the ST in predicting particle properties, the same behavior as ob- 
served at LR = 0.2, which could be attributed to the same physical reasonings discussed in connection 
with the results of that case. 

By comparing Figure 6.1, 1-9 with Figure 6.1. 1-2, two effects of the particles on the gas mean axial velocity 
can be observed. First, the mean axial velocity profile at the inlet plane is flatter than the single-phase 
profile. This is due to the inlet profiles, which correspond to a fully developed pipe flow; the injection 
tube diameter-to-length ratio was equal to 65.2. At this condition, the particles are leading the fluid near 
the injector wall region and thus transfer momentum to the gas, which causes the flattening effect. Sec- 
ond, the mean axial gas velocity downstream of the injector is higher than the single-phase value. For in- 
stance, at z/D = 12.45, an increase of about 20% of the inlet single-phase velocity is caused by the pres- 
ence of the particles. This phenomenon can be attributed to two effects (Mostafa and Elghobashi, 1985). 
One effect is the momentum transfer from the particles to the air, since particle velocity becomes greater 
than gas velocity after a short distance downstream from the injector. The other effect is the modulation 
of the gas turbulence caused by the particles. 

To see how the particles modulate the turbulence structure, the distributions of turbulence kinetic energy 
and shear stress that are shown in Figure 6.1.1-10 can be compared with the corresponding quantities of 
single-phase values illustrated in Figure 6.1. 1-3. At z/D = 12.45, local turbulence intensity is reduced 
from 20% to 12% at the centerline, which corresponds to a reduction of about 40% of the single-phase 
value. This turbulence modulation is caused mainly by the fluctuating relative velocity between the par- 
ticles and the carrier phase. Particles generally cause a reduction in the gas turbulence and an increase in 
the dissipation rate of that energy. This turbulence attenuation reaches its maximum value at a certain 
mass loading ratio, when the particle relaxation time becomes very large compared with the gas La- 
grangian time scale. The performance of the turbulence model for two-phase flows, which considers this 
physical phenomenon, is very good compared with the data in Figures 6.1.1-10 and 6.1.1-11. 

Figure 6.1.1-12 corresponds to measurements of the gas kinetic energy and shear stress at the two load- 
ings and shows the predictions with the single-phase k-e model and its version for two-phase flows 
(Mostafa and Mongia, 1987). In the latter, the turbulence modulation is simulated by introducing extra 
terms in the turbulence kinetic energy and its dissipation rate equations. Figure 6.1.1-12 shows that the 
single-phase model does not predict the turbulence modulation caused by the particles in the two cases. 
However, the two-phase flow model yields fairly good agreement with the data. This result confirms 
previous findings that the interaction between the gas and particles is indeed due to both relative mean 
and fluctuating motion between the two phases, and in some cases, the turbulence modulation caused by 
the particles becomes equally important to the particle dispersion due to gas turbulence. 
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Figure 6.1. 1-1. Test facility: unconfined injector only. 
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Figure 6.1 .1-2. Radial profiles of normalized mean axial velocity, turbulence kinetic energy, and shear 

stress for the single-phase jet. 
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Figure 6.1. 1-3. Radial profiles of normalized rms velocity components and turbulence kinetic energy for 

the single-phase jet 
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figure 6.1.1-4. Radial profiles of normalized mean axial velocity, turbulence kinetic energy, and shear 

stress for the single-phase jet; influence of initial e. 
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Figure 6.1.1-6. Radial profiles of normalized gas axial velocity and particle axial and radial velocity 

components, LR - 0.2. 
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Figure 6.1. 1-7. Radial profiles of normalized gas turbulence kinetic energy, shear stress, and axial 
turbulence intensity, and partide axial turbulence intensity at loading ratio, LR ** 0.2. 
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Figure 6. 1.1 -8. Radial profiles of rms radial and azimuthal velocity components of gas and particles at 

loading ratio, LR » 0.2. 


416 


mill! 'll II! II 1 MUD III 







Stochastic 

treatment 

Deterministic 

treatment 



"'"c TE92-862 

Figure 6.1 .1-9. Radial profile* of normalized gas axial velocity and particle axial and radial velocity 

components, LR ■ 1.0. 
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Figure 6.1.1-10. Radial profiles of normalized gas turbulence kinetic energy, shear stress, and axial 
turbulence intensity, and particle axial turbulence intensity at loading ratio, LR * 1.0. 
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Figure 6.1.1-11. Radial profiles of rms radial and azimuthal velocity components of gas and particles at 

loading ratio, LR = 1.0. 
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Figure 6.1.1-12. Effect of single-phase and two-phase turbulence models on gas phase turbulence kinetic 
energy and shear stress for the two loading ratios, LR * 0.2 and LR « 1 .0. 
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6T2 Confined Single. Round bt 

In this section, a numerical study of a confined flow condition is presented. The purpose of this study is 
to assess the performance of the second-moment closure for confined axisymmetric single-phase and par- 
tide-laden flows. The spray model is based on the stochastic-Lagrangian approach. A higher order nu- 
merical scheme is employed in the continuous phase. The computational results are compared with de- 
tailed experimental data obtained using a two-component phase Doppler technique. A sketch of the test 
section is shown in Figure 6 . 1 . 2 - 1 . The flow conditions used for the cases considered here are given in 
Table 6.1-n. 

The set of governing partial differential equations applied for nonswirling continuous (carrier) phase flow 
consists of equations for continuity, axial (x) and radial (r) momenta, e, and four Reynolds stress compo- 
nents. The finite volume approach (Patankar, 1980) is used to reduce the continuous equations to a set of 
coupled discrete equations. The numerical solutions are obtained using the flux-spline differendng 
scheme, designated FSDS (Varejao, 1979). In this scheme, the total flux is assumed to vary linearly within 
a computational cell (control volume). This assumption leads to a scheme in which the discretization co- 
efficients are identical to those from the exponential scheme (Patankar, 1980) but there is an additional 
source which involves the differences in fluxes at adjacent faces of a cell. The presence of this source term 
enables the flux-spline scheme to respond to the presence of sources and/or multidimensionality of the 
flow. 

The elliptic nature of transfer equations requires that boundary conditions be specified on the four sides 
of the solution domain. Four kinds of boundaries need consideration: inlet, axis of symmetry, wall, and 
the outlet. At the inner boundary, which is located at the first measurement plane (x = 4 mm), the mea- 
sured profiles of U, V, W, u^, v^, w^, and uv are applied. The inlet dissipation rate is prescribed based on 
the assumption of constant length scale and the turbulent kinetic energy, namely 

e - JC±_ (63) 

0.2 R 3 

where R 3 = 76 mm is the pipe radius. The macrolength scale of 0.2 R 3 was estimated through the sensi- 
tivity analysis. 


Table 6 . 1 -II. 

Confined single ict experimental flow conditions. 


Continuous phase 


Medium 

air 

Density, r 

1.17 (Kg/m 3 ) 

Inner jet mass flow rate, mi 

0.0021 (Kg/s) 

Dilute jet mass flow rate, m 2 

0.0272 (kg/s) 

Averaged velocity of inner jet, U 0 
‘phase 

3.935 (m/s) 

Medium 

glass beads 

Density, rp 

2500 (kg/m 3 ) 

Beads diameter. Dp 

105 (mm) 

Centerline velocity 

4.2 (m/s) 

Centerline beads rale 

90(l/s) 

Mass flow rate 

0.00193 (kg/s) 

Mass loading beads to air, LR 

0.925 
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At the axis of symmetry, the radial velocity, shear stresses, and radial gradients of other variables are set 
to zero. All streamwide gradients are presumed zero in the exit plane of the calculation domain, except 
for axial velocity, which needs to be corrected to satisfy integral mass balance. The conventional loga- 
rithmic law, which is based on the local equilibrium assumption, is applied to determine the wall shear 
stress component. The shear stress is then used as boundary condition for U momentum and uv equa- 
tions. The boundary conditions for normal stresses are imposed through the modifications of their pro- 
duction terms using the new calculated wall shear stress value. 

The coupled equations and boundary conditions are solved numerically in a sequential manner using the 
staggered grids for velocities and shear stresses. The main advantage of staggering the location of 
stresses is to enhance the numerical stability, a result of high coupling between the shear stresses and re- 
lated mean strains. The iteration sequence employs the SIMPLER algorithm (Patankar, 1980) to handle 
the coupling between the continuity and momentum equations. The algebraic equations are solved using 
a line-by-line tridiagonal matrix algorithm (TDMA). 

The ordinary differential equations describing the dispersed phase are solved using forward numerical 
integration. The calculations start from the first measurement plane (x = 4 mm) where the mean and 
fluctuation velocity profiles of particles and their rates of injection are available. At the wall, particles are 
assumed to bounce back with an angle of reflection equal to the angle of incidence. 

Numerical iterations are performed over the continuous and dispersed phases until the absolute sums of 
the normalized mass and momentum of the carrier phase at all internal grid points and the change of the 
particle source terms are less than 10~ 5 . 

The predicted mean and turbulence quantities obtained for a single-phase confined round jet are first 
compared with the experimental data. The results presented here were obtained using a 61 x 57 (x - r) 
grid (Figure 6.1 .2-2). The x and r grid coordinates are also shown in Table 6.1-III. A finer grid spacing is 
assigned near the inlet, centerline, and in shear layer. The computational domain extends from the first 
experimental location (x = 4 mm) to 450 mm downstream of the jet exit. Since the measured flow does 
not show any x dependence at x >200 mm, the specified condition d/dx = 0 at the exit plane of the calcula- 
tion domain is realistic. Since the computational results did not change using a finer grid, it is therefore 
assumed that the predicted results are grid independent. 
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Table 6.1-III. 

Confined single jet experimental flow co nditions. 


STREAMWISE COORDINATES OF THE GRID 


I DX 


X XU 


1 

0 . 000E+00 

0 . OOOE+OO 

0. 000E+00 

2 

7.750E-04 

7.750E-04 

0. OOOE+OO 

3 

1 .625E-03 

2.400E-03 

1 .550E-03 

4 

2.050E-03 

4 . 450E-03 

3 . 250E-03 

5 

2.500E-03 

6.950E-03 

5 . 650E-03 

6 

2 . 700E-03 

9.650E-03 

8 . 250E-03 

7 

2.938E-03 

1 . 259E-02 

1 . 1 05E-02 

8 

3 . 313E-03 

1 . 590E-02 

1.413E-02 

9 

3.888E-03 

1 . 979E-02 

1 . 768E-02 

10 

4 . 588E-03 

2 . 438E-02 

2.190E-02 

11 

5.250E-03 

2 . 963E-02 

2.685E-02 

12 

5.775E-03 

3 . 540E-02 

3.240E-02 

13 

6 . 238E-03 

4. 164E-02 

3 . 840E-02 

14 

6.713E-03 

4.835E-02 

4.488E-02 

15 

7.063E-03 

5.541E-02 

5.183E-02 

16 

7.200E-03 

6 . 26 IE-02 

5.900E-02 

17 

7.250E-03 

6.986E-02 

6.623E-02 

18 

7.288E-03 

7.715E-02 

7 . 350E-02 

19 

7.300E-03 

8 . 445E-02 

8 . 080E-02 

20 

7 . 300E-03 

9 . 1 75E-02 

8.810E-02 

21 

7 . 288E-03 

9 . 904E-02 

9.540E-02 

22 

7.275E-03 

1 . 063E-0 1 

1 .027E-01 

23 

7 . 288E-03 

1.136E-01 

1 .100E-01 

24 

7 . 300E-03 

1 .209E-01 

1 .173E-01 

25 

7 . 300E-03 

1 .262E-01 

1 .246E-01 

26 

7.300E-03 

1 . 355E-01 

1 . 319E-01 

27 

7 . 300E-03 

1 . 428E-0 1 

1 . 392E-01 

28 

7.300E-03 

1 . 50 1E-0 1 

1 .465E-01 

29 

7 . 300E-03 

1 .574E-01 

1 . 538E-01 

30 

7.300E-03 

1 .647E-01 

1 .61 IE-01 

31 

7 . 300E-03 

1 .720E-01 

1 . 684E-01 

32 

7 . 300E-03 

1 .793E-01 

1 . 757E-01 

33 

7 . 300E-03 

1 . 866E-0 1 

1 .830E-01 

34 

7 . 300E-03 

1 .939E-01 

1 .903E-01 

35 

7.300E-03 

2.012E-01 

1 .976E-01 

36 

7 . 300E-03 

2 . 085E-0 1 

2.049E-01 

37 

7.300E-03 

2.158E-01 

2. 1 22E-01 

38 

7.300E-03 

2.231 E-0 1 

2.195E-01 

39 

7.200E-03 

2 . 303E-01 

2.268E-01 

40 

7.200E-03 

2.375E-01 

2.339E-01 

41 

7 . 400E-03 

2.449E-01 

2 . 412E-01 

42 

7.400E-03 

2.523E-01 

2.487E-01 

43 

7.325E-03 

2 . 596E-01 

2.560E-01 

44 

7 . 325E-03 

2.670E-01 

2.633E-01 

45 

7.275E-03 

2.742E-01 

2 . 706E-01 

46 

7.275E-03 

2.815E-01 

2 . 779E-01 

47 

7 . 300E-03 

2.888E-01 

2.852E-01 

48 

7.450E-03 

2.963E-01 

2 . 925E-01 
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Table 6.1 -III (cont). 


49 

7.775E-03 

3. 040E-01 

3.001E-01 


50 

7.975E-03 

3. 120E-01 

3 . O80E-OT 

- 

51 

8 . 500E-03 

3.205E-01 

3.160E-01 


52 

9.500E-03 

3.300E-01 

3,250E-01 

i 

53 

1 .000E-02 

3.400E-01 

3.350E-01 

= 

54 

1 . 000E-02 

3.500E-01 

3.450E-01 


55 

1 . 125E-02 

3.613E-01 

3,550E-01 


56 

1 .375E-02 

3. 750E-01 

3 675E-0T 


57 

1 .500E-02 

3 . 900E-01 

825E-01 

- 

58 

1 .500E-02 

4.050E-01 

3. 975E-0T 

- 

59 

1 . 500E-02 

4.200E-01 

4 . 125E-01 


60 

1 .875E-02 

4 . 388E-01 

4.275E-01 

- 

61 

1 . 125E-02 

4.500E-01 

4.500E-01 



TRANSVERS 

COORDINATES OF 

THE GRID 


J 

DY 

Y 

YV 

- 

1 

0 . OOOE+OO 

0. 000E+00 

0. 000E+00 


2 

3 . 750ET-04 

3.750E-04 

0. 000E+00 

= 

3 

7.500E-04 

1 . 125E-03 

7.500E-04 


4 

8 . 750E-04 

2. 000E-03 

1 .500E-03 


5 

1 .000E-03 

3. 000E-03 

2.500E-03 

- 

6 

1 .000E-03 

4 ".'000E-03 

3.500E-03 


7 

1 .000E-03 

5. OOOE-03 

4 * 500E-03 


8 

1 .OOOE-D3 

6 . 000E-03 

5.500E-03 


9 

1 .OOOE-03 

7. 000E-03 

6 . 500E-03 

— 

10 

T , OOOE-03 

8. 000E-03 

7.500E-03 


11 

1 , OOOE-03 

9. OOOE-03 

8.500E-03 


12 

1 .000E-03 

1 - 000E-02 

9 . 500E-03 


13 

1 . 000E-03 

1 . 100E-02 

1 .050E-02 

-- 

14 

1 . 000E-03 

1 .200E-02 

1 .150E-02 


15 

1 .OOOE-03 

1 .300E-02 

1 .250E-02 


16 

1 .OOOE-03 

1 .400E-02 

1 . 350E-02 


17 

1 .000E-03 

1 . 500E-02 

1 .450E-02 


18 

1 . 000E-Q3 

1 . 600E-02 

1 . 550E-02 


19 

1 . 000E-03 

1 .700E-02 

1 . 650E-02 


20 

1 .000E-03 

1 *800E*02 

1 .750E-02 


21 

1 .000E-03 

1 . 900E-G2 

1 . 850E-02 


22 

1 . OOOE-03 

2.000E-02 

1 . 950E-02 


23 

1 .000E-03 

2. 100E-02 

2. 050E-02 


24 

1 .OOOE-03 

2 . 200E-02 

2 . 150E-02 


25 

1 . OOOE-03 

2.300E-02 

2.250E-02 


26 

1 . 000E-03 

2.4OOE-02 

2 . 350E-02 


27 

1 , 000E-03 

2.500E-02 

2.450E-02 


28 

1 . 000E-03 

2 . 600E-02 

2.550E-02 


29 

1 .OOOE-03 

2 . 700E-02 

2.650E-02 


30 

1 .000E-03 

2 . 800E-02 

2.750E-02 


31 

1 .OQOE-03 

2.900E-02 

2.850E-02 


32 

1 . OOOE-03 

3 . 000E-02 

2.950E-02 
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Tabic 6.1 -III (cont). 


33 

1 .000E-03 

3. 100E-02 

3.050E-02 

34 

1 . OOOE-03 

3 .200E-02 

3.150E-02 

35 

1 .250E-03 

3.325E-02 

3.250E-02 

36 

1 .750E-03 

3.500E-02 

3.400E-02 

37 

2.000E-03 

3.700E-02 

3.600E-02 

38 

2.000E-03 

3.900E-02 

3 . 800E-02 

39 

2.000E-03 

4.100E-02 

4.000E-02 

40 

2.000E-03 

4.300E-02 

4 . 200E-02 

41 

2.000E-03 

4.500E-02 

4.400E-02 

42 

2.000E-03 

4.700E-02 

4.600E-02 

43 

2 . 000E-03 

4.900E-02 

4.800E-02 

44 

2.000E-03 

5. 100E-02 

5 . 000E-02 

45 

2.000E-03 

5.300E-02 

5 .200E-02 

46 

2.000E-03 

5 . 500E-02 

5.400E-02 

47 

2.000E-03 

5.700E-02 

5.600E-02 

48 

2.000E-03 

5.900E-02 

5.800E-02 

49 

2.000E-03 

6. 100E-02 

6.000E-02 

50 

2.000E-03 

6.300E-02 

6.200E-02 

51 

2.000E-03 

6.500E-02 

6.400E-02 

52 

2.000E-03 

6 . 700E-02 

6.600E-02 

53 

2.000E-03 

6.900E-02 

6.800E-02 

54 

2 . 000E-03 

7.100E-02 

7.000E-02 

55 

2.000E-03 

7.300E-02 

7.200E-02 

56 

2.000E-03 

7 . 500E-02 

7.400E-02 

57 

1 .000E-03 

7.600E-02 

7.600E-02 
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The present calculations have been made using the flux-spline differencing scheme and the calculated re- 
sults are essentially free of numerical diffusion. Therefore, the discrepancy between the experimental 
data and the prediction can be attributed to two sources, improper boundary conditions at the inlet plane 
and the deficiency of the turbulence model. As regards the inlet conditions, all quantities except the dis- 
sipation rate, e, were prescribed from the measurement. These profiles are shown in Figure 6.1. 2-3 and 
have been normalized by the inner jet averaged velocity, U 0 , and the pipe radius, R3. Above r/R3 > 0.5, 
the flow field is similar to the plug flow condition created by the strong coflow. This region is not the fo- 
cus of this study and is not shown in the presented figures. The uncertainties in the derivation of inlet £ 
profile would adversely affect the calculation at downstream locations. It was shown, however, that the 
use of inlet e derived from the constant length-scale assumption can result in a better prediction (Nikjooy 
and Mongia, 1991). 

A comparison of the normalized axial mean velocity profiles at various locations with data is shown in 
Figure 6.1.2-4. These results were obtained from the k-e model. The predictions arc in good agreement 
with measurements. The velocity profiles are flat near the centerline. The experimental data indicate that 
the centerline velocity constantly decreases. The deceleration in axial velocity is due to the pressure ef- 
fects. Comparisons of the turbulent shear stress and kinetic energy are presented in Figures 6.1 .2-5 and 
6.1.2-6. The k-e model overestimated the maximum level of the turbulent shear stress near the nozzle 
exit, however, the differences decreased in downstream locations. The agreement between the predicted 
and experimental values of kinetic energy is not as good as that for the axial velocity. Even though the 
trends are similar, the predicted kinetic energies are smaller than those derived from the measurements in 
downstream regions. The mean velocity profiles predicted by k-e and DSM are also very close, except for 
some minor differences near the centerline (Figure 6.1 .2-7). 
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Comparisons of the mean and turbulence quantities predicted by the algebraic second-moment (ASM) 
closure with the measured values are presented in Figures 6.1.2-8 through 6.1.2-11. Effects of the pres- 
sure-strain correlation model on mean and turbulence flow fields are analyzed. Two different model con- 
stants for Ci and C2 have been examined in this study. These are Cl = 3.0, C2 = 0.3 (Gibson and Younis, 
1986) and Cl = 5.0, C2 = 0.0 (Donaldson, 1969). The prediction of mean axial velocity profiles using these 
two models are very similar. There are substantial differences between the measured profiles and the 
ASM calculation which show more mixing, and hence lower velocity, near the centerline. 

The predicted stress components are shown in Figures 6.1 .2-9, 6.1.2-10, and 6.1.2-11. Despite the similari- 
ties that appeared in the mean velocity field, the Reynolds stress profiles are different in the shear layer 
and near the centerline zone. The discrepancy observed in the turbulence field between the prediction 
and the measurement is related in large to the deficiency of the pressure-strain model. At downstream 
locations, the differences between the calculations have vanished. The turbulence intensities become al- 
most uniform near the centerline indicating a well-mixed flow condition. 

Comparisons of the normalized mean axial velocity and Reynolds stresses predicted by DSM and ASM 
closures with the measured values are presented in Figures 6. T.2-12 through 6.1.2-16. The pressure redis- 
tribution model used for these calculations involved only the retum-to-isotropy part (4>ij2 = 0) with a con- 
stant Cl = 5.0 in the model. The calculated velocity field predicted by both models is very similar at 
all axial locations. Some minor differences are observed between the two models at downstream regions. 
Overall, the prediction is in good agreement with data. Examination of the calculated axial rms profiles 
indicates that the DSM's results mimic the experimental data better as flow proceeds towards down- 
stream regions. The maximum radial and tangential rms predicted by both models are very close and are 
slightly underpredicted despite the differences that appeared in their axial components. As regards the 
turbulent shear stress, the calculated profiles are similar to the exhibited data trend, however, the k-e 
model resulted in better prediction (Figure 6.1.2-18). 

The calculation for a particle-laden jet was also performed over the same computational domain. The 
same grid distribution was used. As regards the inlet conditions, all quantities except the dissipation rate, 
e, were prescribed from the measurement. These profiles are shown in Figures 6.1.2-19 and 6.1.2-20 and 
have been normalized by the inner jet averaged velocity, U G , particles flux at the centerline, Nq, and the 
pipe radius, R3. The inlet dissipation rate is prescribed based on the assumption of constant length scale 
and the turbulent kinetic energy (Equation 63). 

Comparisons of normalized mean axial velocity, turbulent shear stress, and kinetic energy profiles pre- 
dicted by k-e and DSM closure with the measured values are presented in Figures 6.1.2-21, 6.1.2-22, and 
6.1.2-23. The predicted velocity result is in excellent agreement with the experimental data. The velocity 
profiles are quite flat near the centerline. The centerline velocity decreases for about 10% from the inlet 
plane to station x = 200 mm. The deceleration in axial velocity probably results from the pressure effect. 
The predicted pressure distribution shows a negative radial gradient favoring an outward motion. As a 
result, the axial velocity is slowing down to satisfy the mass conservation. The predicted turbulent shear 
stress and kinetic energy profiles are somewhat larger than those calculated by DSM closure. The k-e 
model overpredicted the maximum uv and k to station x = 35 mm, however, DSM underpredicted the 
peaks at all the streamwise stations. Overall, the turbulence quantities were qualitatively well predicted 
and their behaviors are in agreement with the measured profiles, despite the differences that appeared in 
the shear layer. 

Comparisons of the mean and turbulence quantities predicted by the ASM and DSM closures with the 
experimental data are shown in Figures 6.1.2-24 through 6.1 .2-28. The calculated mean field obtained 
from both models is essentially identical and is in excellent agreement with data. The differences be- 
tween the prediction and data in the downstream region (x > 150 mm) could be related to the diffusion 
process. This is cleared by noting the underprediction of Reynolds stresses in those regions. Compar- 
isons of the predicted and measured streamwise turbulence intensity profiles show some minor differ- 
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ences between the ASM and DSM in calculating the peak values. The differences appear to be more sig- 
nificant for single-phase flow than this two-phase case. At most of the locations the turbulence intensity 
near the centerline was underpredicted. This is due in part to an excessive turbulence dissipation rate at 
the inlet. Nevertheless, the tuibulence intensity profiles obtained from these models are similar to the ex- 
hibited experimental trend. The discrepancy observed in the turbulence field between the prediction and 
measurement is related in large to the deficiency of the pressure-strain model. This has resulted in un- 
derprediction of the maximum values of radial and tangential turbulence intensities. The discrepancies 
between the prediction and measurement are less significant in the case of shear stresses. In the down- 
stream region (x > 150 mm), the calculated results deviate from the data in the outer zone. The results 
demonstrate that the dispersed phase had no significant effects on the continuous phase. 

Figure 6.1.2-29 presents the mean velocity of the dispersed phase. In this calculation the two-phase k-e 
turbulence model was used. It can be seen that the ST approach provides good predictions compared 
with the experimental data, whereas the DT treatments perform poorly for the particles' velocity. Accord- 
ing to the DT approach, a particle moves radially due to its initial mean radial velocity and/or the mean 
radial gas velocity, both of which are very small compared with the axial component. This might explain 
the narrow distribution of particle mean axial velocity predicted by the DT. 

Predicted and measured mean particle velocities predicted by the DSM and ST approach along the axis 
are shown in Figure 6.1.2-30. The calculated velocities arc based on the particles' mass flow rate weighted 
averages. Only a portion of the entire radial section where particles could be found are presented. The 
predicted axial values are in good agreement with data. Despite the variations observed in the gas phase 
velocity near the centerline, the particle velocity shows almost no change. This is related to the fact that 
the particle dynamic relaxation time, xj, is very large compared to the turbulent characteristic time, x. 

The reason for large id is found in the large particle-to-gas density ratio (2500:1.17). This clearly shows 
that the inertia force is responsible for the particles' motion and the effects of drag force are marginal. The 
predicted profiles of mean radial velocity of particles are in qualitative agreement with the data. The ra- 
dial velocity component, however, is at least one order of magnitude smaller than the axial component. It 
is, therefore, reasonable to believe that their quantitative disagreement would not severely affect other re- 
sults. 

The agreements between the predicted and measured fluctuating velocities of the particles are not favor- 
able (Figures 6.1.2-31, 6.1.2-32, and 6.1.2-33). The data show similar radial and tangential fluctuating 
components. On the other hand, the axial rms value is about 1.5 times larger than others. The model has 
successfully predicted the anisotropy feature of the dispersed flow field. However, all three components 
are predicted to be higher than the data. 
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Figure 6.1 2 - 2 . Confined round jet grid layout. 
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4 mm 
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I turbulent shear stress by k-e model with data. 
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Figure 6.12-6. Comparison of predicted turbulent kinetic energy by k-e model with data. 







axial velocity by Jc-c and D9M with data. 
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Figure 6.1 2-1 1 . Comparison of calculated tangential turbulence intensity by ASM with data. 
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Figure 6.1.2-13. Comparison of calculated streamwise turbulence intensity by DSM and ASM with data. 
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Figure 6.1.2-18. Comparison of calculated turbulent shear stress profiles by k-e, ASM, and DSM with 
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Figure 6.12-20. Measured profiles 
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figure 6.12-22. Comparison of predicted gas phase shear stress profiles by DSM and k-e with data, 

particle-laden jet 
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Figure 6 . 1 . 2 - 24 . Comparison of predicted gas phase U velocity profiles by ASM and DSM with data, 

particle-laden jet 
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Figure 6.1.2-25. Companson of calculated gas phase streamwise turbulence intensity profiles by ASM 

and DSM with data, particle- laden ijet. 
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Figure 6.12-27. Comparison of calculated gas phase tangential turbulence intensity profiles by ASM and 

DSM with data, particle-laden jet. 












Figure 6.12-32. Comparison of calculated dispersed phase radial turbulence intensity by DSM with data. 
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6.2 UNCONFINED SINGLE ANNULAR ]ET 


In this section, study of an unconfined annular flow configuration (Figure 6.2-1 ) is presented. The flow 
conditions used for the case considered here are given in Table 6.2-1. In this configuration, the test section 
was directed vertically downward and the annular jet discharged into stagnant air. Measurements have 
been taken using a two-color, two-component laser anemometer system at eight axial stations: 3, 15, 25, 
35, 50, 75, 150, and 300 mm from the exit plane of the jet exit. At each spatial point, the laser simultane- 
ously measured two orthogonal components of velocity. To get all three components, two scans were 
taken. One was used to measure U, V, u 2 , v 2 , and uv components and the other was able to measure U, 
W, u 2 , w 2 , and uw components. Thus all three mean and rms velocity components, uv and uw were 
measured along with U, and u 2 was measured twice. 

The computational mesh used for all calculations consisted of 61 x 70 nonuniformly distributed grid 
points in the axial (x) and radial (r) directions (Figure 6.2-2). The tabulated x and r grid coordinates are 
shown in Table 6.2-II. A finer grid spacing was used near the inlet, centerline, and in the shear layer. The 
computational domain extended from the first measurement plane to 450 mm downstream of the nozzle 
exit. In the radial direction, the entrainment boundary was placed at 170 mm from the axis of symmetry. 
The convergence criterion used to terminate the iterations was the absolute sums of the mass and momen- 
tum residuals at all internal grid points, normalized by inlet mass and momentum fluxes, less than 10 . 

A calc ula tion procedure for elliptic flow requires boundary conditions on all boundaries of the computa- 
tional domain. Four kinds of boundaries need consideration, namely, inlet, axis of symmetry, outlet, and 
the entrainment boundary. At the inlet boundary, which was located at the first measurement plane (x = 
3 mm), the measured profiles of U, V, u 2 , v 2 , w 2 , and uv were applied. These profiles are shown in Fig- 
ure 6.2-3. Two different approaches were used to prescribe the inlet dissipation rate. The first approach 
is based on the assumptions of constant length scale and turbulent kinetic energy, namely 

e = k 1 - 5 /0.20Dj (64) 

where Dj is the jet diameter. The second approach, based on the assumption of local equilibrium (Pk = 
pe) for turbulence energy, can be written as 


dr 

At the axis of symmetry, the radial velocity, shear stresses, and radial gradients of other variables are set 
to zero. At the outlet, axial diffusion is neglected for all variables. Along the entrainment boundary, 
which was placed sufficiently far from the axis of symmetry, the quantity (rV) was assumed zero, and k 
and e were assigned arbitrary low values yielding an eddy viscosity, pt = lOp. 


Table 6.2-1. 

Unconfined single annular jet experimental flow con ditions. 


Ri 

R2 

h 

Uo 

Annular flow rate 


13.82 (mm) 
18.35 (mm) 
7.0 (mm) 

6.9 (m/s) 
0.0033 (Kg/s) 
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Table 6.2-II. 

Single annular jet grid definition. 


STREAMWISE COORDINATES OF THE GRID 


I 

DX 

X 

XU 

1 

0 . OOOE+OO 

0 . OOOE + OO 

o.oooE+oa 

2 

7.750E-04 

7.750E-04 

0 . 000E+00 

3 

1 . 625E-03 

2.400E-03 

1 .550E-03 

4 

2.050E-03 

4.450E-03 

3.250E-03 

5 

2 . 500E-03 

6 . 950E-03 

5.650E-03 

6 

2.700E-03 

9.650E-03 

8 . 250E-03 

7 

2.938E-03 

1 .259E-02 

T.T05E-02 

8 

3.313E-03 

1 .590E-02 

1 .413E-02 

9 

3 . 888E-03 

1 .979E-02 

1 .768E-02 

10 

4.588E-03 

2.438E-02 

2.190E-02 

11 

5.250E-03 

2-963E-02 

2.685E-02 

12 

5.775E-03 

3.540E-02 

3.2A0E-02 

13 

6 . 238E-03 

4 . 164E-02 

3.840E-02 

14 

6.713E-03 

4.835E-02 

4.488E-02 

15 

7.063E-03 

5.541E-02 

5 . 183E-02 

16 

7 . 200E-03 

6.261E-02 

5.900E-02 

17 

7.250E-03 

6.986E-02 

6.623E-02 

18 

7.288E-03 

7.715E-02 

7.350E-02 

19 

7 . 300E-03 

8.445E-02 

8.080E-02 

20 

7.300E-03 

9. 175E-02 

8.810E-02 

21 

7 . 288E-03 

9.904E-02 

9 . 540E-O2 

22 

7 . 275E-03 

1 .063E-01 

1 .027E-01 

23 

7.288E-03 

1 . 136E-01 

1 .100E-01 

24 

7.300E-03 

1 .209E-01 

1 * 173E-01 

25 

7 , 300E-03 

1 . 282E-0 1 

1 .246E-01 

26 

7 . 300E-03 

1 .355E-01 

1 .319E-01 

27 

7 . 300E-03 

1 .428E-01 

1 .392E-01 

28 

7.300E-03 

1 .501E-01 

1 .465E-01 

29 

7.300E-03 

1 .574E-01 

1 . 538E-01 

30 

7.300E-03 

1 . 647E-0 1 

1 .61 IE-01 

31 

7.300E-03 

1 .720E-01 

1 .684E-01 

32 

7, 300E-03 

1 .793E-01 

1 . 757E-01 

33 

7.300E-03 

1 . 866E-0 1 

1 .830E-01 

34 

7.300E-03 

1 .939E-01 

1 .903E-01 

35 

7.300E-03 

2.012E-01 

1 .976E-01 

36 

7.300E-03 

2.085E-01 

2 . 049E-01 

37 

7 . 300E-03 

2. 158E-01 

2.122E-01 

38 

7.300E-03 

2.231E-01 

2 . 195E-01 

39 

7 . 200E-03 

2.303E-01 

2.268E-01 

40 

7 . 200E-03 

2.375E-01 

2.339E-01 

41 

7.400E-03 

2.449E-01 

2.412E-01 

42 

7.400E-03 

2 . 523E-0 1 

2.487E-01 

43 

7.325E-03 

2.596E-01 

2.560E-01 

44 

7.325E-03 

2.670E-01 

2.633E-01 

45 

7.275E-03 

2.742E-01 

2.706E-01 

46 

7.275E-03 

2.815E-01 

2 . 779E-0 1 

47 

7 .300E-03 

2.888E-01 

2 . 852E-01 


PARTI OF 3 
TC82-1140 


462 


Table 6.2-II (cont). 


48 

7.450E-03 

2.963E-01 

2.925E-01 

49 

7.775E-03 

3.040E-01 

3.001E-01 

50 

7.975E-03 

3. 120E-01 

3.080E-01 

51 

8.500E-03 

3 . 205E-01 

3.160E-01 

52 

9 . 500E-03 

3 . 300E-0 1 

3.250E-01 

53 

1 .000E-02 

3.400E-01 

3 .350E-01 

54 

1 .000E-02 

3.500E-01 

3.450E-01 

55 

1 . 125E-02 

3 . 61 3E~0 1 

3.550E-01 

56 

1 . 375E-02 

3.750E-01 

3.675E-01 

57 

1 . 500E-02 

3.900E-01 

3.825E-01 

58 

1 .500E-02 

4 . 050E-0 1 

3 . 975E-01 

59 

1 . 500E-02 

4.200E-01 

4.125E-01 

60 

1 .500E-02 

4.350E-01 

4.275E-01 

61 

7.500E-03 

4 . 425E-01 

4.425E-01 



TRANSVERS 

COORDINATES OF 

THE GRID 

J 

DY 

Y 

YV 

1 

0 . 000E+00 

5.000E-04 

0 . 000E+00 

2 

2 . 500E-04 

7 . 500E-04 

5 . 000E-04 

3 

7.500E-04 

1 . 500E-03 

1 . 000E-03 

4 

1 . 000E-03 

2 . 500E-03 

2 . 000E-03 

5 

1 . 000E-03 

3 . 500E-03 

3 . OOOE-03 

6 

1 . 000E-03 

4.500E-03 

4.000E-03 

7 

1 . 000E-03 

5.500E-03 

5. 000E-03 

8 

1 .000E-03 

6.500E-03 

6. 000E-03 

9 

1 , 000E-03 

7.500E-03 

7. 0O0E-O3 

10 

1 . 000E-03 

8.500E-03 

8. 000E-03 

n 

1 . 000E-03 

9 . 500E-03 

9. 000E-03 

12 

1 . 000E-03 

1 .050E-02 

1 . 000E-02 

13 

1 .000E-03 

1 . 150E-02 

1 . 1 00E-02 

14 

1 .000E-03 

1 .250E-02 

1 . 200E-02 

15 

7 , 500E-04 

1 .325E-02 

1 . 300E-02 

16 

5.000E-04 

1 .375E-02 

1 . 350E-02 

17 

5 , 000E-04 

1 . 425E-02 

1 .400E-02 

16 

5.000E-04 

1 .475E-02 

1 .450E-02 

19 

5.000E-04 

1 . 525E-02 

1 . 500E-02 

20 

5 . 000E-04 

1 . 575E-02 

1 . 550E-02 

21 

5.000E-04 

1 .625E-02 

1 .600E-02 

22 

5 . 000E-04 

1 .675E-02 

1 .650E-02 

23 

5.000E-04 

1 . 725E-02 

1 .700E-02 

24 

5 . O0OE-O4 

1 .775E-02 

1 .750E-02 

25 

5 . 000E-04 

1 .825E-02 

1 .800E-02 

26 

5.000E-04 

1 . 875E-02 

1 .850E-02 

27 

5.000E-04 

1 . 925E-02 

1 . 900E-02 

28 

5.000E-04 

1 . 975E-02 

1 .950E-02 

29 

7.500E-04 

2 . 050E-02 

2.000E-02 

30 

1 . 000E-Q3 

2. 150E-02 

2 . 1 00E-02 

31 

1 .000E-03 

2.250E-02 

2.200E-02 
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Table 6.2-II (eont). 


32 

1 . 00OE-O3 

2.350E-02 

2.300E-02 


33 

1 . 000E-Q3 

2 * 450E-02 

2 . 300E-Q2 


33 

1 .000E-03 

2 , 550E-02 

2.500E-02 

: 

35 

1 . OO0E-03 

2.650E-02 

2.600E-02 

- 

36 

1 .000E-03 

2.750E-02 

2 . 700E-02 


37 

1 .000E-03 

2 . 850E-02 

2.800E-02 


38 

1 . 000E-03 

2 * 950E-02 

2 . 900E-02 

= 

39 

1 .500E-03 

3. 100E-02 

3 . 000E-02 


30 

2.000E-03 

3.300E-02 

3.200E-02 


31 

2.500E-03 

3 . 550E-02 

3.300E-02 


32 

3.000E-03 

3.850E-02 

3.700E-02 


33 

3.000E-03 

3.150E-02 

3 . 000E-02 


33 

3 . 000E-03 

4.450E-02 

3.300E-02 


35 

3.500E-03 

3.800E-02 

4.600E-02 


36 

3 .000E-03 

5 . 200E-02 

5 . 000E-02 

r 

37 

3.500E-03 

5.650E-02 

5.300E-02 

" 

38 

5.000E-03 

6. 150E-02 

5.900E-02 

- 

39 

5 . 000E-03 

6.650E-02 

6 . 300E-02 


50 

5 . 000E-03 

7. 150E-02 

6 . 900E-02 

= 

51 

5.000E-03 

7 . 650E-02 

7 . 300E-02 


52 

5 . OOOE-03 

8. 150E-02 

7 . 900E-02 


53 

5 . 000E-03 

8.650E-02 

8.300E”02 


53 

5. OOOE-03 

9.150E-02 

8 . 900E-02 

= 

55 

5. OOOE-03 

9 . 650E-02 

9.300E-02 


56 

5.500E-03 

1 .020E-01 

9.900E-02 


57 

5.500E-03 

1 .075E-01 

1 .050E-01 


58 

5. 000E-03 

1 . 125E-01 

1 .100E-01 


59 

5. 000E-03 

1 ,175E-01 

1 . 150E-01 


60 

5 . 000E-03 

1 .225E-01 

1 .200E-01 


61 

5 . 000E-03 

1 . 275E-0 1 

1 .250E-01 

_ 

62 

5 - 000E-03 

1 . 325E-0 1 

1 . 300E-01 


63 

5. 000E-03 

1 .375E-01 

1 . 350E-01 



63 

5. OOOE-03 

1 . 325E-0 1 

1 .400E-01 


65 

5. 000E-03 

1 . 375E-0 1 

1 .450E-01 

= 

66 

5. 000E-03 

1 . 525E-0 1 

1 . 500E-01 


67 

5, 000E-03 

1 . 575E-0 1 

1 .550E-01 


68 

5. 000E-03 

1 . 625E-0 1 

1 .600E-01 


69 

5. 000E-03 

1 .675E-01 

1 .650E-01 


70 

2 . 500E-03 

1 .700E-01 

1 .700E-01 



PART 3 OF 3 
TE321140 


In an attempt to assess the importance of inlet turbulence dissipation rate (ej n ) on mean velocity and 
stress fields predictions, two different approaches were followed to determine ej n (Equations 64 and 65). 
Comparisons of the calculated mean axial velocity profiles by DSM with the experimental data are illus- 
trated in Figure 6.2-4. The results show differences in the prediction of peak and centerline values in the 
developing region. The inlet e profile calculated, based on equilibrium assumption (Equation 65), re- 
sulted in a closer prediction with data near the inlet plane. Further downstream, the differences between 
the model and data become more severe. In developing flow region, the experiment shows that the mo- 
mentum is being transferred towards the centerline, increasing the centerline velocity. However, the 
model did not simulate this process adequately. As a result, the velocity was underpredicted in the inner 
region, however, the DSM performed reasonably in the outer region. At downstream (x > 150 mm), the 
calculated velocity field is in good agreement with measurement, although some minor discrepancies ap- 
pear near the centerline. 

Comparisons of the predicted turbulence intensities from the DSM closure are shown in Figures 6.2-5, 6.2- 
6, 6.2-7, and 6.2-8. When the inlet e distribution is calculated with the constant length scale assumption. 


464 


the profiles are closer to the exhibited data trend. A wide disparity exists between this calculation and 
the calculation with ejn obtained from Equation 65. The model indicates a faster mixing in the inner re- 
gion than is shown by experiment. Numerical study indicates that the inlet e profile is a very important 
factor in predicting the maximum values of mean and turbulence quantities. The turbulent normal 
stresses are decreased due to an excessive dissipation rate at the inlet. Comparisons of radial and tangen- 
tial Reynolds stress components reveal that essentially the same trends are observed for v^ and w^ pro- 
files. The effects of the ej n on the shear stress profiles are found to be serious near the inlet region where 
the central recirculation is formed (Figure 6.2-8). In comparison with the data, the calculated values are 
not in good agreement in the developing flow region. A wide disparity exists between the calculations. 
The €j n derived from constant length scale assumption resulted in more realistic profiles. However, the 
locations of the maxima and minima have been shifted in the developing region of the flow field. 

The predicted mean axial velocity profiles from the k-e model are compared with data in Figure 6.2-9. 
These results were also obtained using the two different inlet dissipation rate profiles. Similar behavior is 
also observed here. The major differences between the two calculations appeared in the region where 
maximum velocity occurs. Both methods are unable to predict the flow in the center region accurately. 
But the discrepancies are removed as flow proceeds towards the downstream region. Figure 6.2-10 shows 
the kinetic energy profiles at various axial locations. Again, the behavior of two sets of computations is in 
closer agreement to each other than to measurement. Most of the differences are seen in the peak and 
centerline regions. The agreement between the predicted and experimental values of kinetic energy is not 
as good as that for the axial velocity. Even though to some extent the trends are similar, the predicted ki- 
netic energy levels are smaller than those derived from the measurements. For the ease of comparison, 
the predicted mean axial velocity and turbulent shear stress profiles from the DSM and k-e models have 
been compared in Figures 6.2-11 through 6.2-14. Depending on the shape of ei n distribution, the discrep- 
ancies between the DSM and k-e model predictions can be severe or minor. This clearly demonstrates the 
importance of inlet profiles, particularly the dissipation rate for turbulence model validation (Sturgess et 
al, 1983). 
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Figure 6.2-1. Unconfined annular flow configuration. 
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Figure 6.2-4. Comparison of predicted mean axial velocity by DSM with data (1 of 4). 
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Figure 6.2-4. Comparison of predicted mean axial velocity by DSM with data (2 of 4). 
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Figure 6.2-4. Comparison of predicted mean axial velocity by DSM with data (3 of 4). 
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Figure 6.2-4. Comparison of predicted mean axial velocity by DSM with data (4 of 4). 
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Figure 6.2-5. Comparison of calculated rms axial velocity profiles by DSM with data (1 of 4). 
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Figure 6.2-5. Comparison of calculated rms axial velocity profiles by DSM with data (2 of 4). 
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o 



Figure 6.2-5. Comparison of calculated rms axial velocity profiles Ity DSM with data (3 of 4). 



Figure 6.2-5. Comparison of calculated rms axial velocity profiles by DSM with data (4 of 4). 
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Figure 6.2-6. Comparison of calculated rms radial velocity profiles by DSM with data (1 of 4). 
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Figure 6.2-6. Comparison of calculated rms radial velocity profiles by DSM with data (2 of 4). 
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Figure 6.2-6. Comparison of calculated rms radial velocity profiles by DSM with data (3 of 4). 
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Figure 6.2-6. Comparison of calculated rms radial velocity profiles by DSM with data (4 of 4). 
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Figure 6.2-7. Comparison of calculated rms tangential velocity profiles by DSM with data (1 of 4). 



Figure 6.2-7. Comparison of calculated rms tangential velocity profiles by DSM with data (2 of 4). 
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Figure 6.2-7. Comparison of calculated rms tangential velocity profiles by DSM with data (3 of 4). 
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Figure 6.2-7. Comparison of calculated rms tangential velocity profiles by DSM with data (4 of 4). 
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Figure 6.2-8. Comparison of calculated turbulent shear stress profiles by DSM with data (3 of 4). 



Figure 6.2-8. Comparison of calculated turbulent shear stress profiles by DSM with data (4 of 4). 
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Figure 62-9. Comparison of predicted mean axial velocity by k-e model with data (1 of 4) 
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Figure 62-9. Comparison of predicted mean axial vdocity by k-e model with data (2 of 4) 




Figure 6.2-9. Comparison of predicted mean axial velocity by k-e model with data (3 of 4 ). 



Figure 6.2-9. Comparison of predicted mean axial vdocity by k-e modd with data (4 of 4). 
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: igure 62-10. Comparison of predicted turbulent kinetic energy by k-e model with data (1 of 4). 
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Figure 6.2-10. Comparison of predicted turbulent kinetic energy by k-e model with data (2 of 4). 
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Figure 6.2-10. Comparison of predicted turbulent kinetic energy by k-e model with data (3 of 4). 
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Figure 6.2-10. Comparison of predicted turbulent kinetic energy by k-e model with data (4 of 4). 
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Figure 6.2-11. Comparison of calculated U velocity by k-e and DSM with data (3 of 4). 


o 



Figure 622-11. Comparison of calculated U velocity by k-e and DSM with data (4 of 4). 
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Figure 6.2-12. Comparison of calculated uv profiles by k-e and DSM with data (1 of 4). 



Figure 6.2-12. Comparison of calculated uv profiles by k-e and DSM with data (2 of 4). 







uy i 1/ 


u V'/ a / TE92-1 180 


Figure 62-13. Comparison of calculated U velocity by k-e and D5M with data (1 of 4). 
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Figure 62-13. Comparison of calculated U velocity by k-e and DSM with data (2 of 4). 






Figure 62-13. Comparison of calculated U velocity by lc-e and DSM with data (3 of 4). 


o 



Figure 6.2-13. Comparison of calculated U velocity by k-e and DSM with data (4 of 4). 




Figure 6.2-14. .Comparison of calculated uv profiles by k-e and DSM with data (1 of 4). 
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Figure 6.2-14. Comparison of calculated uv profiles by k-e and DSM with data (2 of 4). 




Figure 6.2-14. Comparison of calculated uv profiles by k-e and DSM with data (3 of 4). 



Figure 6.2-14. Comparison of calculated uv profiles by k-e and DSM with data (4 of 4). 
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III , 1,1, 



63 SINGLE ANNULAR SWIRLING JET 


63,1 Uncon fined Annular Swirling let 

In this section, computations for an unconfined axisymmetric annular swirling jet are reported. This con- 
figuration is the same as the previous one, but with a swirler in the annular passage. In this configura- 
tion, an annular jet was directed vertically downward within a 457 mm^ wire mesh screen. Experimental 
data were obtained at eight axial locations 3, 15, 25, 35, 50, 75, 150, and 300 mm from the exit plane of the 
injector. The flow conditions used for this case are given in Table 63-1 and the sketch of the test section is 
shown in Figure 63.1-1. 

The governing equations were discretized using the control volume approach (Patankar, 1980). The flux- 
spline differencing scheme was applied for obtaining the numerical scheme. The computational mesh 
used for all calculations consisted of 61 x 70 nonuniformly distributed grid points in the axial (x) and ra- 
dial (r) directions. A finer grid spacing was used near the inlet, centerline, and in the shear layer (Figure 
63.1-2). The tabulated grid points in the axial and radial directions are shown in Table 6.3-II. The com- 
putational domain extended from the first measurement plane, located downstream of the nozzle exit at a 
distance of 3.0 mm to 450 mm. 

A calculation procedure for elliptic flow needs boundary conditions on all boundaries of the computa- 
tional domain. Four kinds of boundaries, namely, inlet, axis of symmetry, outlet, and the entrainment 
boundary, need consideration. At the inlet boundary located at the first measurement plane, the mea- 
sured profiles of mean velocities and Reynolds stresses were prescribed. These profiles are shown in Fig- 
ure 63.1-3. The k-profile was obtained from the measured Reynolds normal stress components. The e 
value at the inlet plane (ei n ) was derived from the kinetic energy profile and constant length-scale as- 
sumption. At the axis of symmetry, shear stresses, radial velocity, and radial gradient of other variables 
are set to zero. At the outlet, axial diffusion is neglected for all variables. Along the entrainment bound- 
ary, the quantity (rV) was assumed constant. In addition, the axial velocity (U) and tangential velocity 
(W) were assumed to be zero. Turbulent kinetic energy and its dissipation rate were assigned arbitrary 
low values yielding an eddy viscosity, m = 10p.. 

Predictions of mean axial and tangential velocities obtained from the k-e and DSM are illustrated in Fig- 
ures 63.1-4 and 63.1-5. A wide disparity exists between the models' predictions of axial velocity espe- 
cially near the centerline. There are substantial differences in the capability of the various models to 
promote or hinder formation of the recirculation zone. With regard to the comparison between mea- 
surements and calculations, the trend predicted by the k-e model seems to be better away from the cen- 
terline. The maximum velocity was underpredicted by both models. The results of the DSM closure 
show less radial diffusion process. As a result, in comparison with data, the peak values were all shifted. 
One of the main reasons for discrepancies between the model and data can be related to the uncertainties 
in treating the open boundaries. Ibis is clear when the models are applied for confined swirling flow. 

The k-e model also resulted in slightly better prediction for tangential velocity. An examination of the 
calculated Reynolds stresses (Figures 63.1-6, 63.1-7, 63.1-8) indicates that the relative performance of the 


Table 6.3-1. 

Unconfined annular swirling iet test configuration. 


Ri 

R2 

h 

Uo 

Swirl angle 
Swirl airflow rate 


13.82 (mm) 
18.35 (mm) 
7.0 (mm) 

6.9 (m/s) 

60 deg 

0.0033 (Kg/s) 
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Table 6.3-11. 

Uncpnfined annular swirling iet grid definition. 


STREAMWISE COORDINATES OF THE GRID 


I DX X XU 


1 

0 . 000E+00 

0 . OO0E+OQ 

0 . 000E+00 

2 

7 . 750E-04 

7.750E-04 

0 .OOOE+OO 

3 

1 . 625E-03 

2 . 400E-03 

1 .550E-03 

4 

2.050E-03 

4.450E-03 

3.250E-03 

5 

2.500E-03 

6.950E-03 

5.650E-03 

6 

2 . 700E-03 

9.650E-03 

8.250E-03 

7 

2.938E-03 

1 .259E-02 

1105E-02 

8 

3.313E-03 

1 .590E-02 

1 . 413E-02 

9 

3 . 888E-03 

1 .979E-02 

1 .768E-02 

10 

4.588E-03 

2.438E-02 

2 . 190E-02 

11 

5.250E-03 

2.963E-02 

2 . 685E-02 

12 

5.775E-03 

3 . 540E-02 

3 . 240E-02 

13 

6 . 238E-03 

4. 164E-02 

3 . 840E-02 

14 

6.7fSf-5S 

4.835E-02 

4.488E-02 

15 

7.063E-03 

5 . 541E-02 

5 . 183E-02 

16 

7.200E-03 

6 .26 IE -02 

5.900E-02 

17 

7.250E-03 

6.986E-02 

6 .623E-02 

18 

7.288E-03 

7.715E-02 

7 . 350E-02 

19 

7.300E-03 

8.445E-02 

8.080E-02 

20 

7.300E-03 

9. 175E-02 

8.810E-02 

21 

7.288E-03 

9.904E-02 

9. 540E-02 

22 

7.275E-03 

1 .063E-01 

1 . 027E-01 

23 

7.288E-03 

1 .T3SE-01 

1 . 100E-01 

24 

7.300E-03 

1 .209E-01 

1 . 173E-01 

25 

7.300E-03 

1 .282E-01 

1 . 246E-01 

26 

7.300E-03 

1 . 355E-01 

1 .319E-01 

27 

7.300E-03 

1 .428E-01 

1.392E-01 

28 

7.300E-03 

1 .501E-01 

1 .465E-01 

29 

7.300E-03 

1 . 574E-01 

1 .538E-Q1 

30 

7.300E-03 

1 .647E-01 

1 .61 IE-01 

31 

7.300E-03 

1 . 720E-0 1 

1 . 684E-01 

32 

7 , 300E-03 

1 . 793E-Q1 

1 . 757E-01 

33 

7.300E-03 

1 .866E-01 

1 .830E-01 

34 

7.300E-03 

1 . 939E-01 

1.903E-01 

35 

7 . 300E-03 

2 . 0 12E-01 

1 .976E-01 

36 

7 . 300E-03 

2.085E-01 

2 . 049E-01 

37 

7.300E-03 

2. 158E-01 

2.122E-01 

38 

7.300E-03 

2.231E-01 

2 . 195E-01 

39 

7.200E-03 

2.303E-01 

2.268E-01 

40 

7.200E-03 

2 . 375E-01 

2.339E-0! 

41 

7 , 400E-03 

2.449E-01 

2.412E-01 

42 

7.400E-03 

2.523E-01 

2.487E-01 

43 

7.325E-03 

2 . 596E-01 

2.560E-01 

44 

7.325E-03 

2 . 670E-0 1 

2.633E-01 

45 

7.275E-03 

2.742E-01 

2.706E-01 

46 

7 , 275E-03 

2.815E-01 

2. 779E-01 

47 

7.300E-03 

2.888E-01 

2.852E-01 
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Tabic 6.3-11 (cont). 


48 

7.450E-03 

2.963E-01 

2.925E-01 

49 

7 . 775E-03 

3.040E-01 

3.001E-01 

50 

7.975E-03 

3. 120E-01 

3 . 080E-01 

51 

8.500E-03 

3.205E-01 

3.160E-01 

52 

9 . 500E-03 

3.300E-01 

3.250E-01 

53 

1 .000E-02 

3 . 400E-Q 1 

3.350E-01 

54 

1 .000E-02 

3.500E-01 

3.450E-01 

55 

1 .125E-02 

3.613E-01 

3.550E-01 

56 

1 .375E-02 

3.750E-01 

3.675E-01 

57 

1 . 500E-02 

3.900E-01 

3.825E-01 

58 

1 . 500E-02 

4 . 050E-01 

3.975E-01 

59 

1 .500E-02 

4.200E-01 

4.125E-01 

60 

1 . 500E-02 

4.350E-01 

4 . 275E-01 

61 

7.500E-03 

4.425E-01 

4.425E-01 


TRANSVERS 

COORDINATES OF 

THE GRID 


J 

DY 

Y 

YV 

1 

0 . 000E+00 

5.000E-04 

0 . 000E+00 

2 

2 « 500E-04 

7.500E-04 

5.000E-04 

3 

7.500E-04 

1 .500E-03 

1 , 000E-03 

4 

1 . 000E-03 

2 . 500E-03 

2.000E-03 

5 

1 .000E-03 

3 . 500E-03 

3.000E-03 

6 

1 . 000E-03 

4.500E-03 

4 . 000E-03 

7 

1 . 000E-03 

5 . 500E-03 

5.000E-03 

8 

1 .000E-03 

6 . 500E-03 

6 . 000E-03 

9 

1 . 000E-03 

7 . 500E-03 

7.000E-03 

10 

1 .000E-03 

8.500E-03 

8 . 000E-03 

11 

1 .000E-03 

9 . 500E-03 

9.000E-03 

12 

1 .000E-03 

1 . 050E-02 

1 . OOOE-02 

13 

1 .000E-03 

1 .150E-02 

1 . 100E-02 

14 

1 . 000E-03 

1 . 250E-02 

1 .200E-02 

15 

7.500E-04 

1 . 325E-02 

1 . 300E-02 

16 

5 . Q00E-04 

1 . 375E-02 

1 . 350E-02 

17 

5.000E-04 

1 . 425E-02 

1 . 400E-02 

18 

5.000E-04 

1 . 475E-02 

1 .450E-02 

19 

5.000E-04 

1 . 525E-02 

1 . 500E-02 

20 

5.000E-04 

1 . 575E-02 

1 . 550E-02 

21 

5 . 000E-04 

1 .625E-02 

1 . 600E-02 

22 

5.000E-04 

1 . 675E-02 

1 .650E-02 

23 

5 . 000E-04 

1 . 725E-02 

1 . 700E-02 

24 

5.000E-04 

1 .775E-02 

1 . 750E-02 

25 

5.000E-04 

1 .825E-02 

1 . 800E-02 

26 

5.000E-04 

1 . 875E-02 

1 . 850E-02 

27 

5.000E-04 

1 . 925E-02 

1 . 900E-02 

28 

5 . 000E-04 

1 . 975E-02 

1 . 950E-02 

29 

7.500E-04 

2 . 050E-02 

2 . 000E-02 

30 

1 .000E-03 

2. 150E-02 

2.100E-02 

31 

1 .000E-03 

2.250E-02 

2.200E-02 
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Table 6.3-II (cont). 


32 

1 .OOOE-03 

2.350E-02 

2.300E-02 

33 

1 . OOOE-03 

2.450E-02 

2.400E-02 

34 

1 . OOOE-03 

2.550E-02 

2.S00E-02 

35 

1 .000E-03 

2.650E-02 

2.600E-02 

36 

1 .000E-03 

2T75OE-02 

2 . 700E-02 

37 

1 . OOOE-03 

2.850E-02 

2.800E-02 

38 

1 .OOOE-03 

2 , 950E-02 

2^950F-02 

39 

T.500E-03 

3. i00E-02 

3.000E-02 

40 

2. OOOE-03 

3.300E-02 

3.200E-02 

41 

2 . 500E-03 

3 . 550E-02 

3.400E-02 

42 

3. 000E-03 

3.850E-02 

3.700E-02 

43 

3 . O00E-O3 

4.150E-02 

4 , 0O0E-02 

44 

3. OOOE-03 

47450E-02 

A.300E-02 

45 

3.500E-03 

4 . 800E-02 

4.600E-02 

46 

4. OOOE-03 

5.200E-02 

5.000E-02 

47 

4 . 500T-03 

5.650E-02 

5.400E-02 

48 

5.000E-03 

6. 150E-02 

5.900E-02 

49 

5.000E-03 

6.650E-02 

6 . 400E-02 

50 

5. OOOE-03 

7.150E-02 

6.900E-02 

51 

5 . OOOE-03 

7.650E-02 

7.400E-02 

52 

5. 000E-03 

8.150E-02 

7.900E-02 

53 

5. OOOE-03 

8.650E-02 

8.400E-02 

54 

5. OOOE-03 

9. 150E-02 

8.900E-02 

55 

5 . OOOE-03 

9 . 650E-02 

9.400E-02 

56 

5 . 500E-03 

1 . 020E-01 

9.900E-02 

57 

5.500E-03 

1 .075E-01 

1 .050E-01 

58 

5. 000E-03 

1 . 125E-01 

1 . 100E-01 

59 

5. 000E-03 

1 . 175E-01 

1 . 150E-01 

60 

5 . 000E-03 

1 .225E-01 

1 .200E-01 

61 

5. 000E-03 

1 .275E-01 

1 .250E-01 

62 

5. 000E-03 

1 . 325E-0 1 

1 .300E-01 

63 

5. 000E-03 

1 .375E-01 

1 .350E-01 

64 

5 . 000E-03 

1 .425E-01 

1 .400E-01 

65 

5 .000E-03 

1 . 475E-01 

1 .450E-01 

66 

5. 000E-03 

1 . 525E-0 1 

1 .500E-01 

67 

5. OOOE-03 

1 . 575E-0 1 

1 . 550E-01 

68 

5. OOOE-03 

1 . 625E-01 

1 .600E-01 

69 

5. 000E-03 

1 .675E-01 

1 .650E-01 

70 

2.500E-03 

1 . 700E-0 1 

1 . 700E-01 
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model is dependent on the flow region. The results show differences, especially in the prediction of the 
centerline values. The discrepancy can be the result of the inlet dissipation rate profile. The inlet dissipa- 
tion rate is a determining factor in predicting the maximum level of turbulence intensity. Comparison of 
the calculated shear stress profiles with data is shown in Figure 6.3.1-9. In comparison with the mea- 
surement, the predicted values by the k-e model are in better agreement, however, the maximum and 
minimum have not been well predicted . 

In an attempt to assess the importance of inlet dissipation rate (ejn) on flow field, calculations were also 
made using an alternative ej n distribution, which was derived from local equilibrium assumption for tur- 
bulence energy (Equation 65). Comparison of the calculated mean axial and tangential velocity profiles 
by the k-e model with the experimental data are presented in Figures6.3.1-10 and 6.3.1-11. Numerical ex- 
periments indicate that the inlet e profile is the single most important factor in predicting the maximum 
values of mean and turbulence quantities. This point is now clear by comparing the two sets of results. 
The major differences between these two conditions are in the peak region. The inlet dissipation rate 
profile derived from the equilibrium assumption resulted in shorter maximum velocity and smaller recir- 
culation zone to satisfy the global mass conservation. The tangential velocity has also been affected simi- 
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larly by ei n . Figures 6.3.1-12 and 6.3.1-13 show comparisons of tuibulent shear stress and kinetic energy 
profiles. The ein found by equilibrium assumptions resulted in higher shear and energy. These will in- 
crease the radial diffusion process, and as a result, the peak velocities will be quickly diminished. 
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TE92-1 189-4 

Figure 6.3.1 -1. Unconfined annular swirling Jet geometrical details. 
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figure 63.1-3. Measured profiles at the inlet (x = 3 mm) (2 of 2). 
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Figure 63.1-4. Comparison of calculated mean axial velocity by DSM and k-e with data (1 of 3). 



Figure 63.1-4. Comparison of calculated mean axial velocity by DSM and k-e with data (2 of 3). 






10.0 20.0 30.0 40.0 50.0 a 0.0 10.0 20.0 30.0 40.0 



gure 63.1-5. Comparison of calculated mean tangential velocity by DSM and k-e with dab * (2 of 3). 



Figure 63.1-5. Comparison of calculated mean tangential velocity by DSM and k-e with data (3 of 3). 
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Figure 63,1-6. Comparison of calculated axial rms velocity by DSM with data (1 of 3). 
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Figure 63.1-6. Comparison of calculated axial rms velocity by DSM with data (2 of 3). 








Figure 63.1-6. Companion of calculated axial rms velocity by DSM with data (3 of 3). 



Figure 63.1-7. Comparison of calculated radial rms velocity by DSM with data (1 of 3). 
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Figure 63.1-7. Comparison of calculated radial rms velocity by D5M with data (2 of 3). 



Figure 63.1-7. Comparison of calculated radial rms velocity by DSM with data (3 of 3). 
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Figure 63.1-8. Comparison of calculated tangential rma velocity by DSM with data (1 of 3). 
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Figure 63.1-8. Comparison of calculated tangential nns velocity by DSM with data (2 of 3). 
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Figure 63.1-8. Comparison o f calculated tangential rms velocity by DSM with data (3 of 3). 



Figure 6.3.1-9. Comparison of calculated turbulent shear stress profiles by k-e and DSM with data (1 of 3) 
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Figure 6.3.1 -9. Comparison of calculated turbulent shear stress profiles by k-e and DSM with data (3 of 3). 




Figure 63.1-10. Comparison of predicted mean axial velocity by k-e with data (1 of 3). 
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Figure 63.1-10. Comparison of predicted mean axial velocity by k-e with data (2 of 3). 







Figure 6.3.1-10. Comparison of predicted mean axial velodly by k-e with data (3 of 3). 



Figure 63.1-1 1 . Comparison of predicted mean tangential velocity by k-e with data (1 of 3). 





Figure 63.1-11. Comparison of predicted mean tangential velocity by k-e with data (2 of 3). 
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Figure 63.1-11. Comparison of predicted mean tangential velocity by k-e with data (3 of 3). 





Figure 6-3-1-12. Comparison of predicted turbulent shear ffress with d#ta(l of 3). 



Figure 6.3.1-12. Comparison of predicted turbulent shew stress with data (2 of 3). 
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Figure 63.1-13. Comparison of predicted turbulent kinetic energy with data (2 of 3). 
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Figure 6.3.1-13. Comparison of predicted turbulent kinetic energy with data (3 of 3). 




I 


The predicted mean and turbulence fields obtained from the k-e, ASM, and DSM closures are compared 
with the experimental data. The experiments were performed in a modeled laboratory combustor and 
the complex flow field has an aerodynamically controlled, swirl stabilized recirculation zone. A two- 
color, two-component laser anemometer system was used to measure the velocity components. At each 
spatial point, the laser simultaneously measured two orthogonal components of velocity. To get all three 
components, two scans were taken. One was used to measure U, V, u^, v^, and uv components, and the 
other one was able to measure U, W, u^, w^, and uv components. Thus all three mean and fluctuating 
velocity components were measured, with U and u^ measured twice. Data were obtained at eight axial 
stations 4, 15, 25, 35, 50, 75, 100, and 150 mm from the exit of the swirler. At each station, between 10 and 
20 radial points were scanned, as determined by the desired level of profile resolution. The experimental 
test parameters are summarized in Table 6.3-III and the sketch of the test section is shown in Figure 6.3.2- 
1 . 


A nonuni form grid of 48 x 44 in the axial and radial directions was used for all calculations (Figure 6.3.2- 
2). A finer grid spacing was used near the inlet, centerline, and in the shear layer. The streamwise and 
transverse coordinates of the grid points are shown in Table 6.3-IV. The elliptic nature of the transfer 
equations requires that boundary conditions be specified on the four sides of the solution domain. Four 
kinds of boundaries need consideration: inlet, axis of symmetry, wall, and outlet. At the inlet boundary, 
which is located at the first measurement plane (x = 4 mm), the measured profiles of U, V, W, u^, v^, w , 
uv, and uw are applied. These profiles are shown in Figure 6.3.2-3. Three different approaches were 
used to prescribe the inlet dissipation rate. The first approach is based on the assumption of constant 
length scale and turbulent kinetic energy, namely 



( 66 ) 


where Dj is the jet diameter. Here, the macrolength scale of 0.2Dj was estimated through the sensitivity 
analysis. The second approach is derived from the k-e eddy viscosity and neglecting the 3V/8X term at 
the inlet to relate the turbulent shear stresses to their corresponding mean strain rates. That is. 


e =Cuk 


/auf /aw wfT 

t [l 3r / [dr r) 
[(uv) 2 +(vw) 2 ] 1/2 


(67) 


Table 6.3-III. 


R1 

R2 

R3 

h 

Uo 

Swirl angle 
Swirl airflow rate 
Dilute airflow rate 


13.82 mm 
1835 mm 

76.00 mm 

7.0 mm 
6.9 m/s 
60 

0.0033 kg/s 
0.272 kg/s 


515 


Table 6.3-IV. 

Confined annular swirling jet grid definition. 


STREAMWISE COORDINATES OF THE GRID 


I DX X XU 


1 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

2 

7.750E-04 

7.750E-04 

O.OOOE+OO 

3 

1 .625E-03 

2 . 400E-03 

1 .550E-03 

4 

2 . 050E-03 

4 . 450E-03 

3.250E-03 

5 

2.500E-03 

6.950E-03 

5.650E-03 

6 

2.700E-03 

9.650E-03 

8 . 250E-03 

7 

2 . 938E-03 

1 . 259E-02 

1 .105E-02 

8 

3.313E-03 

1 .590E-02 

1 . 413E-02 

9 

3.888E-03 

1 . 979E-02 

1 .768E-02 

10 

4 . 588E-03 

2.438E-02 

2.190E-02 

1 1 

5.250E-03 

2 * 963E-02 

2.685E-02 

12 

5 . 775E-03 

3.540E-02 

3.240E-02 

T3 

6 . 238E-03 

4. 164E-02 

3.840E-02 

14 

6.713E-03 

4 . 835E-02 

4 . 488E-02 

15 

7 . 063E-03 

5.541E-02 

5.183E-02 

16 

7 . 200E-03 

6.261E-02 

5.900E-02 

17 

7 . 250E-Q3 

6.9S6E“0r 

0623E-P2 

18 

7.288E-03 

7.715E-02 

7 . 350E-02 

19 

7 . 300E-03 

8.445E-02 

8.080E-02 

20 

7 , 300E-03 

9. 175E-02 

8.810E-02 

21 

7 , 288E-03 

9.904E-02 

9.540E-02 

22 

7 . 275E-03 

1 .063E-01 

1 .027E-01 

23 

7 . 288E-03 

1 . 136E-01 

1.100E-01 

24 

7 , 300E-03 

1 .209E-01 

1 .173E-01 

25 

7 . 30GE-03 

1 .282E-01 

1 .246E-01 

26 

7.300E-03 

1 .355E-01 

1 .319E-01 

27 

7.300E-03 

1 .428E-01 

1 .392E-01 

28 

7.300E-03 

1 .501E-01 

1 .465E-01 

29 

7.300E-03 

1 . 574E-0 1 

1 .538E-01 

30 

7.300E-03 

1 . 647E-0 1 

1 .611E-01 

31 

7.300E-03 

1 .720E-01 

1 .684E-01 

32 

7.300E-03 

1 -793E-01 

1 . 757E-0 1 

33 

7.300E-03 

1 .866E-01 

1 .830E-01 

34 

7.300E-03 

1 .939E-01 

1 . 903E-01 

35 

7 . 300E-03 

2.012E-01 

1 .976E-01 

36 

7.300E-03 

2.085E-01 

2.049E-01 

37 

7.300E-03 

2 . 1 58E-01 

2 . 122E-01 

38 

7.300E-03 

2.231E-01 

2.195E-01 

39 

7 . 200E-03 

2.303E-01 

2.268E-01 

40 

7 . 200E-03 

2.375E-01 

2.339E-01 

41 

7.400E-03 

2.449E-01 

2.412E-01 

42 

7.400E-03 

2.523E-01 

2.487E-01 

43 

7.325E-03 

2.596E-01 

2.560E-01 

44 

7.325E-03 

2 . 670E-0 1 

2.633E-01 

45 

7 . 275E-03 

2.742E-01 

2.706E-01 

46 

7 . 275E-03 

2.815E-01 

2.779E-01 

47 

7 . 300E-03 

2.888E-01 

2.852E-01 

48 

3 . 650E-03 

2.925E-01 

2.925E-01 
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Table 6.3-IV (cont). 

TRANSVERS COORDINATES OF THE GRID 


J 

DY 

Y 

YV 

1 

0 . OOOE+OO 

0 , 000E+00 

0. 000E+00 

2 

5 . 000E-04 

5 . OOOE-04 

0. 000E+00 

3 

1 .500E-03 

2. 000E-03 

1 .000E-03 

4 

2 . 000E-03 

4. 000E-03 

3 . OOOE-03 

5 

2.000E-03 

6. 000E-03 

5. 000E-03 

6 

2 . 000E-03 

8. 000E-03 

7. 000E-03 

7 

2.000E-03 

1 .000E-02 

9. 000E-03 

8 

2.000E-03 

1 . 200E-02 

1 .100E-02 

9 

2 . 000E-03 

1 . 400E-02 

1 . 300E-02 

10 

1 .750E-03 

1 .575E-02 

1 .500E-02 

11 

1 .250E-03 

1 . 700E-02 

1 .650E-02 

12 

1 .000E-03 

1 , 800E-02 

1 . 750E-02 

13 

1 .000E-03 

1 . 900E-02 

1 . 850E-02 

14 

1 . 000E-03 

2 . 000E-02 

1 .950E-02 

15 

1 . 000E-03 

2 . 1 OOE-02 

2 . 050E-02 

16 

1 . 250E-03 

2.225E-02 

2 . 150E-02 

17 

1 . 750E-03 

2 . 400E-02 

2.300E-02 

18 

2.000E-03 

2.600E-02 

2 . 500E-02 

19 

2.000E-03 

2 . 800E-02 

2.700E-02 

20 

2 . OOOE-03 

3.000E-02 

2.900E-02 

21 

2.000E-03 

3 . 200E-02 

3 . 1 00E-02 

22 

2. OOOE-03 

3 . 400E-02 

3 . 300E-02 

23 

2. OOOE-03 

3 . 600E-02 

3 . 500E-02 

24 

2. 000E-03 

3 . 800E-02 

3.700E-02 

25 

2. 000E-03 

4 . 000E-02 

3 . 900E-02 

26 

2. 000E-03 

4.200E-02 

4.100E-02 

27 

2. 000E-03 

4.400E-02 

4.300E-02 

28 

2. 000E-03 

4 . 600E-02 

4 . 500E-02 

29 

2. 000E-03 

4 . 800E-02 

4 . 700E-02 

30 

2 , 000E-03 

5 . 000E-02 

4 . 900E-02 

31 

2. 000E-03 

5 . 200E-02 

5 . 1 00E-02 

32 

2. 000E-03 

5 . 400E-02 

5.300E-02 

33 

2. 000E-03 

5 . 600E-02 

5 . 500E-02 

34 

2. 000E-03 

5 . 800E-02 

5 . 700E-02 

35 

2. 000E-03 

6 . OOOE-02 

5 . 900E-02 

36 

2. 000E-03 

6 . 200E-02 

6 . 1 00E-02 

37 

2. 000E-03 

6.400E-02 

6 . 300E-02 

38 

2. 000E-03 

6 . 600E-02 

6 . 500E-02 

39 

2. 000E-03 

6 . 800E-02 

6 . 700E-02 

40 

2. OOOE-03 

7.000E-02 

6 . 900E-02 

41 

2. 000E-03 

7.200E-02 

7.100E-02 

42 

2. OOOE-03 

7.400E-02 

7. 300E-02 

43 

1 . 500E-03 

7.550E-02 

7 . 500E-02 

44 

5.000E-04 

7.600E-02 

7 . 600E-02 
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where C^, = 0.09 is a model constant. The third approach, based on the assumption of local equilibrium 
for turbulence energy (Pk = pe), Boussinesque approximation, isotropic effective viscosity, and ignoring 
stream wise velocity gradient, can be written as 



(68) 
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At the axis of symmetry, the radial velocity, shear stresses, and radial gradients of other variables are set 
to zero. All streamwise gradients are presumed zero in the exit plane of the calculation domain except for 
axial velocity, which needs to be corrected to satisfy integral mass balance. The conventional logarithmic 
law, which is based on the local equilibrium assumption, is applied to determine the wall shear stress 
components in the axial and tangential directions. These shear stresses are used as boundary conditions 
for U and W momentum, as well as for the uv and rw equations. The boundary conditions for normal 
stresses are imposed through the modification of their production terms using the new calculated wall 
shear stresses. 

The computational domain extends from the first experimen tal location (x = 4 mm) to 300 mm down- 
stream of the jet exit. The reason for this selection was to prevent recirculation at the exit Since the mea- 
sured flow does not show any x dependence at x > 150 mm, the specified condition, d/dx = 0, at the exit 
plane of the calculation domain is realistic. The coupled equations and boundary conditions are solved 
numerically in a sequential manner using the staggered grids for velocities and shear stresses (Figure 
63.2-4). The main advantage of staggering the locations of stresses is to enhance numerical stability, a re- 
sult of high coupling between the shear stresses and related mean strains. 

The present calculations have been made using the higher order numerical scheme, and the calculated re- 
sults are essentially free of numerical diffusion. Therefore, the discrepancies between the experimental 
data and the predictions can be attributed to two sources: improper boundary conditions at the inlet 
plane and the deficiencies of the turbulence model. As regards the inlet conditions, all quantities except 
the dissipation rate (e) were prescribed from the measurement. Uncertainties in the derivation of the e 
profile, however, would adversely affect the calculations at downstream locations. In an attempt to as- 
sess the importance of inlet dissipation rate (ej n ) on flow field prediction, three different approaches 
(Equations 66, 67, and 68) were followed to determine e values from the measured quantities. The inlet e 
profiles for these cases are shown in Figure 63.2-5 and have been normalized by the inlet mean velocity 
(Uo) and pipe radius (R3), 

Overall, the major differences among the inlet dissipation profiles appeared in the peak and the cento-line 
regions. Comparisons of the predicted mean velocities and turbulence quantities obtained from the k-e 
model with measurement are shown in Figures 63.2-6 through 63.2-9. It is noted that the use of inlet 
dissipation calculated from the strain rate and shear stress resulted in higher peak and centerline veloci- 
ties in upstream planes. However, the recirculation region has been better predicted as flow proceeds 
towards the downstream region. There are significant differences among the calculated turbulent kinetic 
energy and shear stress profiles. The main discrepancy has arisen in prediction of the maximum values 
in die shear layer. 

Comparisons of the predicted mean velocities and Reynolds stresses profiles using the DSM closure are 
presented in Figures 6.3.2-10 through 6.3.2-13. These profiles have been normalized by the inlet velocity 
and pipe radius. Above r/R3 > 0.5, the flow field is similar to plug flow created by the strong coflow. 

This region is not the focus of our study and is not shown in the figures. The use of ejn profiles derived 
from the constant length-scale assumption (Equation 66) and the equilibrium assumption (Equation 68) 
resulted in similar predictions outside the recirculation zone. However, the calculated flow field with ej n 
obtained from the shear stress and strain rate (Equation 67) showed somewhat different behavior. This is 
related mainly to the maximum value of the inlet e, which has been changed considerably. The variations 
of the ejn profiles near the centerline also affect the prediction of the reverse flow region, although the 
Reynolds stresses and dissipation rate are significantly smaller near the centerline, compared to the shear 
region. The uncertainties in the evaluation of velocity gradients may lead to errors in ej n values as well. 

It must be pointed out that, in a similar nonswirling jet configuration, the flow field was found to be less 
sensitive to the magnitude of ej n . The differences between the ejn calculated by various methods were 
larger for the swirling flow case than the nonswirling case. The effects of ejn on shear stress are serious 
near the inlet, where the central recirculation is formed. In comparison with the data, the locations of the 
maximum and minimum have been shifted in the developing region of the flow field. In the outer part of 
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the recirculation zone, where the effect of coflow is more pronounced, and near the axis of symmetry, 
where the velocity gradient is negligible in the radial direction, the discrepancies among the calculations 
are in fact diminished. The peak values are the main differences that appear using various inlet e distri- 
butions. When the inlet e are found from Equations 66 and 68, the predicted turbulence intensity profiles 
are similar to the exhibited data trend. However, a wide disparity exists between these calculations and 
the calculations with ein obtained from Equation 67. The turbulent normal stresses are decreased due to 
an excessive dissipation rate at the inlet. This results in an insufficient diffusion effect and therefore an 
inaccurate prediction of the velocity profile. 

Effects of the pressure-strain correlation model on mean and turbulence fields are also analyzed. The in- 
let dissipation rate is calculated using a constant length-scale assumption (Equation 66). The relative 
merits of the four different pressure redistribution models on ASM closures have been examined in detail 
by Nikjooy and So (1987). Their results indicate that the inability of the second-moment closure to 
correctly resolve the turbulence field could be attributed, to a large extent, to the modeling of the <bij term. 
For swirling flows, the redistribution of turbulence by fluctuating pressure is bound to be different from 
that for the corresponding nonswirling flows. A vast range of values have been suggested for Cl and C2 
by different researchers. However, most choices of Ci and C2 gather close to the line (1 - C2VCI = 0.23 
(Launder, 1989). For a simple shear flow in local equilibrium condition the Reynolds stresses are a func- 
tion of (1 - C2)/C1 only. The original work of Naot et al (1970) neglected the return to isotropy process 
(<t>iji model constant). Conversely, Donaldson's proposal (1969) totally ignores the $ij2 process and just 
includes a contribution from the process, with a large model constant of Ci = 5. Recently, Gibson and 
Younis (1986) argued that, for swirling flows, the constants for 4>ij models should be modified. They pro- 
posed a simple modification of model constants, Ci = 3.0 and C2 = 03, a choice that emphasizes <|>jji more 
than This choice is based on the observation that swirling flows are better predicted with this com- 
bination. These values were used for the calculations shown in Figures 6.3.2-10 through 6.3.2-13. 

Three different model constants for Ci and C2 have been examined in this study. These are Ci = 3.0 and 
C2 = 0.3 (Gison and Younis, 1986), Ci = 1.8 and C2 = 0.6 (Launder, 1989), and Q = 5.0 and C2 = 0.0 
(Donaldson, 1969). The predictions of mean axial and tangential velocity profiles using these three mod- 
els are shown in Figures 633-14 and 6.3.2-15. There are substantial differences between the measured 
profiles and the calculations with the second model, which shows less mixing near the centerline, causing 
the reverse flow to persist. There is, correspondingly, a slowing down in the diminution of the swirl ve- 
locity. As a result, unrealistic recirculation length has been predicted. The mean velocity fields predicted 
by the Gibson and Younis' model and Donaldson's model are almost identical and are in complete accord 
with experiment near the centerline, while the swirl velocity is well predicted at most of the locations. 
With a constant C2 - 0 (<fcj2 = 0), the results obtained for Ci = 1.5 - 10 were all nearly the same. But with 
C2 * 0, convergence was obtained only for certain combinations of Ci and C2 satisfying (1 - C2VC1 = 0.2 
- 0.35. However, the results were different. This seems to suggest that for C2 = 0 results are insensitive to 
Ci, but when C2 * 0 the combination of Ci and C2 is important. It must be pointed out that the slow part 
of the pressure-strain rate term by itself did not fully account for all of the stress-component interactions. 
Therefore selection of an appropriate model for 44j is crucial for swirling flow calculations. A more so- 
phisticated version of <|>jj2 is the linear quasi-isotropic (QI) model (Equation 22). This model (Equation 22) 
includes both the symmetric and antisymmetric mean strain effects on redistribution modeling. The QI 
model achieves reasonable success in thin shear flows but shows no significant advantage over the 
isotropization production (IP) model to swirling flows. Comparison of mean axial velocity profiles that 
resulted from the IP and QI models at location of 35 mm is shown in Figure 6.3.2-16. The solutions were 
obtained for Cl — 1.8 and C2 = 0.6. Some minor improvements are observed near the centerline. How- 
ever, the QI model still gives inferior results in comparison with the Gibson and Younis' model. This in- 
dicates that the additional linear terms, in the present form, did not result in superior predictions. Similar 
observations were made for tangential velocity and Reynolds stress components. 
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The predicted stress components are shown in Figures 63.2-17 through 6.3.2-20. Despite the significant 
differences that appeared in the mean velocity fields, the Reynolds stress profiles are nearly the same. 
This can be explained by noting that the momentum equations consist of the Reynolds stress gradients. 
Although the stresses look similar, careful examination reveals differences in their gradients. Thus any 
variation in the Reynolds stress gradients will result in significant changes in mean velocities. To 
verily this point, a calculation was made using the Reynolds stresses obtained horn one solution (Ci = 5 
and C 2 = 0) and the mean velocities from another solution (Cl = 1.8 and C2 = 0.6). By freezing the stress- 
es' field and solving only the momentum equations, the velocity field corresponding to Cl = 5.0 and C2 = 
0 was obtained. This clearly shows that the differences in the Reynolds stress gradients are responsible 
for differences in the velocity field. 

Comparisons of the mean and turbulence quantities predicted by the ASM and DSM closures with the 
measured values are presented in Figures 6.3.2-21 through 63.2-24. The results Have been shown for four 
stream wise locations, namely, x = 15, 25, 35, and 50 mm. The pressure redistribution model used for 
these calculations involved only the retum-to-isotropy part (foj2 = 0) with a constant Ci * 5.0 in the fon 
model. In general, the results of the ASM and DSM are very similar in the near field (x ~15 mm). At sub- 
sequent downstream locations, there are substantial differences between the two models in the prediction 
of centerline velocity and the recirculation length. Use of the ASM promotes earlier decay of the center- 
line velocity relative to the DSM and fails to display the size of the experimental recirculation zone. The 
reason is attributed to the turbulent diffusion process, which depends on the Reynolds stress gradients. 
Despite the fact that the predicted normal stresses look similar, their gradient profiles are different. In the 
outer flow region, the ASM and DSM predict similar axial velocity distributions. With regard to the 
comparison between measurements and calculations, the predictions by ASM seem to be slightly better 
than DSM for tangential velocity. Examination of the calculated u^ profiles indicates that the ASM's re- 
sults mimic the experimental data better as flow proceeds toward the downstream region. The maximum 
turbulence intensity predicted by both models is very close and is slightly underpredicted by both models 
only at x = 25 mm. As regards the turbulent shear stresses, uv and uw, the calculated profiles are similar 
to the exhibited data trend. The measurement shows a constant increase of the uw values near the center- 
line from x = 15 mm to x = 75 mm, except at x = 25 mm. This sudden change is considered to be an exper- 
imental error. Overall, the shear stresses were qualitatively well predicted and their behaviors are in 
agreement with the measured profiles, although the maximum values were shifted slightly. In general, 
apart from the centerline velocity profile, the ASM closure performed a reasonable job in predictions of 
the tangential velocity and Reynolds stress components. 

The present study serves to illustrate three very important points in the development of turbulence clo- 
sures for calculation of complex swirling flows. First, the inlet e profile was found to be highly influential 
on the stress field prediction. The mean velocity field, to some extent, is affected by the choice of tin es- 
pecially in the recirculation zone. Second, the model assumptions for the pressure-strain correlations are 
not very satisfactory and need improvement. The retum-to-isotropy part of the ftj by itself could predict 
important features of the flow. However, the rapid part, foj2, was found to be influential for swirling 
flow, and the emphasis should be placed on developing a better combination of foji and $jj2- Finally, the 
ASM closure appeared to be unsuitable for prediction of the central recirculation zone. The ASM applica- 
tion, however, resulted in good tangential velocity and stresses profiles. In general, for such cases, a DSM 
closure is more suitable. 
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Figure 6.3.2-3. Measured profiles at the inlet plane (x ■ 4mm). 
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Ftgiire 6.3.2-4. Staggered arrangement of Reynolds stress components. 



Figure 6.3.2-5. Inlet dissipation profiles. 
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Figure 63.2-6. Comparison of calculated mean axial velocity by k-e model using various ej n (1 of 4). 
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Figure 6.32-6. Comparison of calculated mean axial velocity by k-e model using various ei n (2 of 4). 
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Figure 632-7. Comparison of calculated mean tangential velocity by k-e model using various e in (1 of 4) 



Figure 63.2-7. Comparison of calculated mean tangential velocity by k-e model using various (2 of 4). 
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Figure 63.2-7. Comparison of calculated mean tangential velocity by k-e model using various ei n (4 of 4). 
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Figure 63.2-8. Comparison of calculated turbulent kinetic energy by k-e model using various ej n (1 of 4). 




Figure 632-8. Comparison of calculated turbulent kinetic energy by k-e model using various ej n (2 of 4). 
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Figure 63.2-8. Comparison of calculated turbulent kinetic energy by k-e model using various ej n (3 of 4). 




Figure 633-8. Comparison of calculated turbulent kinetic energy by k-e model using various ej n (4 of 4). 
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63.2-9. Comparison of calculated turbulent shear stress by k-e model using various ej n (1 of 4) 
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Figure 63.2-9. Comparison of calculated turbulent shear stress by k-€ model using various z in (2 of 4) 
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Figure 632 - 9 . Comparison of calculated turbulent shear stress by k-e model using various ej n (4 of 4). 
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Figure 63.2-11, Comparison of calculated mean tangential velocity by DSM using various 










Figure 63*2-14. Comparison of calculated mean axial velocity using various models with data. 



Figure 632-15. Comparison of calculated mean tangential velocity using various <>, models with data. 
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Figure 6.3.2-16. Comparison of calculated U velocity with IP and QI models for confined annular 

swirling jet. 
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Figure 63.2-19. Comparison of calculated radial turbulence intensity with 








Figure 63.2-20. Comparison of the predicted mean and turbulence quantities 
by ASM and DSM with data. 
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Figure 6.3.2-21. Comparison of the predicted mean and turbulence quantities 
by ASM and DSM with data. 
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Figure 6.3.2-23. Comparison of the predicted mean and turbulence quantities 
by ASM and DSM with data. 




6.4 COAXIAL JETS 

644 Unwnfined Costal fete 


For the present study, an unconfined flow configuration was selected, operating with an axial jet injector 
surrounded by a nonswirling annular jet, as shown in Figure 6.4.1-1. In this configuration, the injector 
was directed vertically downward within a 457 mm^ wire mesh screen. The flow conditions used for the 
case considered are given in Table 6.4-1. The jet center (diameter, D = 24.1 mm) is surrounded by an an- 
nular jet with the inner and outer diameters of 29 mm and 36.7 mm, respectively. The effective area ratio 
and axial velocity ratio ofthe annular jet to the center jet are 0.87 and 1.8, respectively. Data were ob- 
tained at eight axial stations 2, 15, 25, 35, 50, 75, 150, and 300 mm from the exit plane of the injector. At 
each axial station, between 10 and 20 radial points were scanned as determined by the desired level of 
profile resolution. 

The computations for coaxial jets can be made using a parabolic marching procedure if the radial pressure 
gradients are small. Such a situation occurs if velocities in the two streams are comparable or the inner 
streams are faster and if the swirl is weak. However, if the swirl is strong and/or the outer stream is sig- 
nificantly faster, the radial pressure gradients become significant and a region of reverse flow develops. 
The ultimate goal is to extend this study to swirling flow analyses. Therefore, a calculation procedure 
based on elliptic flows was used. The discretization equations are obtained using a control volume ap- 
proach (Pa tankar, 1 980). Three different numerical schemes wer e tested for the convection and diffusion 
terms in the transport equations. The schemes used are the hybrid scheme (Patankar, 1980), the power- 
law differencing scheme (Patankar, 1980), and the flux-spline differencing scheme (Varejao, 1979). 

A calculation procedure for elliptic flow requires boundary conditions on all boundaries of the computa- 
tional domain. Four kinds of boundaries need consideration, namely, inlet, axis of symmetry, outlet, and 
the entrainment boundary. At the inlet boundary, which was located at the first measurement plane, the 
measured profiles of U and V were prescribed. The k-profile was obtained from the measured Reynolds 
stresses. TTiese profiles are shown in Figure 6.4. 1-2. This kinetic energy distribution and the measured 
shear stress profile were used to derive the e values at the inlet plane through the following relationship 


c _Cjik 2 pU\ 
(- u v) \ dr I 


(69) 


Table 6.4-1. 

Unconfined coaxial jets experimental conditions at 2 mm downstream of pipe exit. 


Parameter 
Inner jet 

Centerline velocity (m/s) 

Density, pi (kg/m 3 ) 

Mass flow rate, mi (kg/s) 

Reynolds number. Re = 4mi/7cmD* 

Ratio of particle-to-gas mass flow rate, LR 
Annular jet 

Maximum velocity (m/s) 

Density, pi (kg/m 3 ) 

Mass flow rate, m 2 (kg/s) 

*D = 0.0241 m 


Single- 

Partide-laden flow 

phase flow 

Air 

Particles 

4.64 

4.76 

4.20 

1.178 

1.178 

2500.0 

0.0021 

0.0021 

0.00042 

6000 

6000 



0.2 


6.60 

6.60 


1.178 

1.178 


0.0033 

0.0033 
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At the axis of symmetry, the radial velocity and the radial gradients of other variables are set to zero. At 
the outlet, axial diffusion is neglected for all variables. Along the entrainment boundary, which was 
placed sufficiently far from the axis of symmetry, the quantity (rV) was assumed constant. In addition, 
the axial velocity U was assumed zero and k and e were assigned arbitrarily low values yielding an eddy 
viscosity, m = 10 p. 

The computational mesh used for all calculations consisted of 76 x 69 nonuniformly distributed gnd 
points in the axial (x) and radial (r) directions (Figure 6.4. 1-3). The tabulated x and r grid coordinates are 
shown in Table 6.4-IL A finer grid spacing was used near the inlet, centerline, and in the shear layer. The 
computational domain extended from the first measurement plane, located downstream of the nozzle exit 
at a distance of 2.0 mm, to 40 inner jet diameters downstream of the nozzle exit. In the radial direction, 
the entrainment boundary was placed at a distance of six jet diameters from the axis of symmetry. The 
convergence criterion used to terminate the iterations was the absolute sums of the mass and momentum 
residuals at all internal grid points, normalized by inlet mass and momentum fluxes, less than 5 x 10“^. 
The numerical results were obtained using two turbulence models with various differencing schemes. 

The calculated mean and turbulence quantities are compared with the measurement at selected stations. 

6.4. 1.1 The k-e Turbulence Model 

The effect of different discretization schemes is shown by comparing the predicted axial velocity profiles 
at three stream wise locations, namely, x = 15, 35, and 75 mm. The velocity profiles are presented in Fig- 
ure 6.4.1-4. It is noted that, except for some minor differences, all three schemes for the convective terms 
yield nearly identical results. In earlier studies (Patankar et al, 1987; Varejao, 1979), it was shown that in 
the region of high Peclet number the flux-spline results are more accurate than those from the power-law 
scheme. The fact that for the present situation there are no significant differences between results from 
these schemes indicates that the results are grid independent. The differences between the hybrid and the 
power-law schemes are attributed to the different treatments of the diffusion terms. The computed re- 
sults at the selected axial stations compare reasonably well with the experimental data. The computations 
consistently show sharper gradients than the experiment at the points of the maximum and minimum 
velocity. 

Figure 6.4.1-5 shows the kinetic energy profiles at three axial locations. The experimental kinetic energy 
profiles were derived from the measured Reynolds stresses. In these figures, results from the power-law 
scheme and the flux-spline scheme have been shown. Again, the two sets of computations are in close 
agreement with each other. Most of the differences are seen in the regions of steep gradients where the 
flux-spline results are expected to be more accurate. The agreement between the predicted and experi- 
mental values of kinetic energy is not as good as that for the axial velocity. Even though the trends are 
similar, the predicted kinetic energy levels are smaller than those derived from the measurements. 

Since the present calculations are essentially free of numerical diffusion, the discrepancies between the 
experimental data and the predictions can be attributed to two sources, improper boundary conditions at 
the inlet plane and the deficiencies of the turbulence model. As regards the inlet conditions, all quantities 
except the dissipation rate were prescribed from the experiment. The e values, however, were derived 
from the measured shear stresses and the mean velocity gradients. The uncertainties in the measure- 
ments and in the evaluation of the velocity gradients may lead to errors in the e values which would ad- 
versely affect the calculations at downstream locations. 

Numerical experiments indicate that the inlet e profile is the single most important factor in predicting 
the maximum values of mean and turbulence quantities, provided a reasonable inlet kinetic energy dis- 
tribution is available. To study the sensitivity to inlet e profile, calculations were also made using an al- 
ternative e distribution, which were derived from the turbulent kinetic energy and an assigned length 
scale distribution (3% of the radius). The inlet e profiles for both cases are shown in Figure 6.4.1-6. The 
major differences between these two conditions are near the centerline region, however, the peak values 
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Table 6.4-II. 

Confined coaxial iets grid definition. 


STREAMWISE COORDINATES OF THE GRID 


I 

DX 

X 

XU 

1 

0 .OOOE+OO 

0 . OOOE+OO 

0 . OOoE+OO 

2 

1 .050E-03 

1 .050E-03 

0. OOOE+OO 

3 

2.200E-03 

3 . 250E-03 

2. 100E-03 

4 

2.4DtnEMJ3 

5 . 650E-03 

4.400E-O3 

5 

2.600E-03 

8 . 250E-03 

6.900E-03 

6 

2.800E-03 

1 . 105E-02 

9.60QE-03 

7 

3.075E-03 

1 . 4 1 2E-02 

1 . 250E-02 

8 

3.550E-03 

1 . 767E-02 

1 .575E-02 

9 

4 .225E-03 

2. 190E-02 

1 .960E-02 

10 

4.950E-03 

2 . 685E-02 

2.420E-02 

1 1 

5.550E-03 

3.240E-02 

2.950E-02 

12 

6 . 000E-03 

3.840E-02 

3.530E-02 

13 

6.475E-03 

4.487E-02 

4. 150E-02 

14 

6.950E-03 

5. 182E-02 

4.825E-02 

15 

7. 175E-03 

5.900E-02 

5.540E-O2 

16 

7.225E-03 

6.622E-02 

6.260E-02 

17 

7.275E-03 

7.350E-02 

6 . 985E-02 

18 

7.300E-03 

8 . 080E-02 

7.715E-02 

19 

7.300E-03 

8.810E-02 

8.445E-02 

20 

7.300E-03 

9 . 540E-02 

9 . 175E-02 

21 

7.300E-03 

1 .027E-01 

9.905E-02 

22 

7.300E-03 

1 . 100E-01 

1 .063E-01 

23 

7 .300E-03 

1 . 173E-01 

1 . 136E-01 

24 

7.300E-03 

1 .246E-01 

1 . 209E-0 1 

25 

7.300E-03 

1 .319E-01 

1 . 282E-0 1 

26 

7.300E-03 

1 . 392E-0 1 

1 . 355E-0 1 

27 

7.300E-03 

1 .465E-01 

1 .428E-01 

28 

7.300E-03 

1 .538E-01 

1 .501E-01 

29 

7.300E-03 

1 .61 IE-01 

1 .574E-01 

30 

7.300E-03 

1 .684E-01 

1 . 647E-0 1 

31 

7.300E-03 

1 .757E-01 

1 . 720E-01 

32 

7.300E-03 

1 .830E-01 

1 .793E-01 

33 

7 . 300E-03 

1 .903E-01 

1 .866E-01 

34 

7.300E-03 

1 . 976E-0 1 

1 . 939E-0 1 

35 

7 , 300E-03 

2.049E-Q1 

2.012E-01 

36 

7.300E-03 

2. 122E-01 

2 . 085E-Q1 

37 

7.300E-03 

2. 195E-01 

2. 158E-01 

38 

7.300E-03 

2 . 268E-0 1 

2.231E-01 

39 

7.300E-03 

2.341E-01 

2.304E-01 

40 

7.300E-03 

2.414E-01 

2.377E-01 

41 

7.300E-03 

2.487E-01 

2 . 450E-0 1 

42 

7.300E-03 

2.560E-01 

2.523E-01 

43 

7.300E-03 

2.633E-01 

2.596E-01 

44 

7.300E-03 

2 . 706E-0 1 

2.670E-01 

45 

7.300E-03 

2.779E-01 

2.742E-01 

46 

7.300E-03 

2.852E-01 

2.815E-01 

47 

7.350E-03 

2.925E-01 

2 . 888E-0 1 

48 

7 . 575E-03 

3.001E-01 

2.962E-01 
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Table 6.4-II (cont). 


49 

7 . 875E-03 

3.080E-01 

3.040E-01 

50 

8 .OOOE-03 

3 . 160E-01 

3. 120E-01 

51 

8.000E-03 

3.240E-01 

3 . 200E-0 1 

52 

8 . 000E-G3 

3.320E-01 

3.280E-01 

53 

8.500E-03 

3.405E-01 

3.360E-01 

54 

9.500E-03 

3 . 500E-0 1 

3.450E-01 

55 

1 .0O0E-O2 

3.600E-01 

3.550E-01 

56 

1 . OOOE-02 

3.700E-01 

3.650E-01 

57 

1 . 250E-02 

3.825E-01 

3 . 750E-0 1 

58 

1 . 750E-02 

4.000E-01 

3.900E-01 

59 

2. 000E-02 

4.200E-01 

4 . 100E-01 

60 

2. 000E-02 

4.400E-01 

4 .300E-01 

61 

2. 000E-02 

4 . 600E-0 1 

4 . 500E-01 

62 

2. 000E-02 

4 . 800E-0 1 

4.700E-01 

63 

2.250E-02 

5.025E-01 

4.900E-01 

64 

2.750E-02 

5 . 300E- 0 1 

5. 150E-01 

65 

3. 000E-02 

5.600E-01 

5.450E-01 

66 

3. 000E-02 

5.900E-01 

5 . 750E-0 1 

67 

3 . 250E-02 

6.225E-01 

6.050E-01 

68 

3.750E-02 

6 . 600E-0 1 

6.400E-01 

69 

4. OOOE-02 

7 . 000E-01 

6 . 800E-0 1 

70 

4. 000E-02 

7 . 400E-0 1 

7.200E-01 

71 

4. OOOE-02 

7 . 800E-0 1 

7 . 600E-0 1 

72 

4 . OOOE-02 

8.200E-01 

8.000E-01 

73 

4 . 000E-02 

8.600E-01 

8.400E-01 

74 

4 . OOOE-02 

9.000E-01 

8.800E-01 

75 

5. OOOE-02 

9.500E-01 

9.200E-01 

76 

3. OOOE-02 

9.800E-01 

9.800E-01 


TRANSVERS 

COORDINATES OR 

THE GRID 


J 

DY 

Y 

YV 

1 

0 .000E+00 

0 . OOOE + OO 

0 .OOOE+OO 

2 

3 . 750E-04 

3 . 750E-04 

0 .OOOE+OO 

3 

7.500E-04 

1 . 125E-03 

7 . 500E-04 

4 

8 . 750E-04 

2.000E-03 

1 .500E-03 

5 

1 . 000E-03 

3 . 000E-03 

2.500E-03 

6 

1 . 000E-03 

4 . 00OE-03 

3.500E-03 

7 

1 .000E-03 

5 . 00OE-03 

4 . 500E-03 

8 

1 . 000E-03 

6.000E-03 

5 . 500E-03 

9 

1 . 0O0E-O3 

7 . Q00E-03 

6 . 500E-03 

10 

1 . 000E-03 

8 . 000E-03 

7.500E-03 

1 1 

1 . 000E-O3 

9.000E-03 

8 . 500E-03 

12 

1 . OOOE-03 

1 . OOOE-02 

9.500E-03 

13 

1 .000E-03 

1 . 100E-02 

1 . 050E-02 

14 

8.750E-04 

1 . 187E-02 

1 . 1 50E-02 

15 

6.250E-04 

1 . 250E-02 

1 .225E-02 

16 

5 . OOOE-04 

1 .300E-02 

1 . 275E-02 


PART 2 OF 3 
TE92-1246 


551 



Table 6.4-II (cont). 


17 

5.000E-04 

1 .350E-02 

1 .325E-02 

18 

6.250E-04 

1 .412E-02 

1 .375E-02 

19 

8.750E-04 

1 .500E-02 

1 . 450E-02 

20 

1 . 00QE-03 

1 .600E-02 

1 .550E-02 

21 

1 . 000E-03 

1 . 700E- 02 

1 .650E-02 

22 

1 . 000E-Q3 

1 . 800E-02 

1 .750E-02 

23 

1 , OOOE-03 

1 . 900E-Q2 

1 . 850E-02 

24 

9.800E-04 

1 .998E-02 

1 . 950E-02 

25 

1 * 047E-03 

2. 103E-02 

2.046E-02 

26 

1 .300E-03 

2 . 233E-02 

2. 1S9E-02 

27 

1 “ 563E-03 

2.389E-02 

2 . 306E-0Z 

28 

1 . 735E-03 

2 . 562E-02 

2.472E-02 

29 

1 . 875E-03 

2.750E-02 

2 . 653E-02 

30 

1 .995E-03 

2 . 950E-Q2 

2.847E-02 

31 

2. 100E-03 

3, 159E-02 

3.052E-02 

32 

2. 192E-03 

3 . 379E-02 

3 , 267E-02 

33 

2 . 275E-03 

3 . 606E-02 

3 . 490E-02 

34 

2.353E-03 

3.841E-02 

3.722E-02 

35 

2 . 425E-03 

4.084E-02 

3.961E-02 

36 

2 . 495E-03 

4.333E-02 

4 . 207E-02 

37 

2 .557E-03 

4 .589E-02 

4 460E-02 

38 

2 .615E-03 

4 .851E-02 

4.719E-02 

39 

2 . 672E-03 

5. 118E-02 

4.983E-02 

40 

2.725E-03 

5.390E-02 

5.253E-02 

41 

2 . 777E-03 

5.668E-02 

5.528E-02 

42 

2 . 830E-03 

5.951E-02 

5 . 808E-02 

43 

2.877E-03 

6 , 239E-02 

6.094E-02 

44 

2 . 922E-03 

6.531E-02 

6.384E-02 

45 

2.967E-03 

6 . 828E-02 

6 . 678E-02 

46 

3.01 0E-03 

7. 129E-02 

6.977E-02 

47 

3.050E-03 

7 . 434E-02 

7 . 280E-02 

48 

3 ♦ 090E-03 

7.743E-02 

7.587E-02 

49 

3. 128E-03 

8.056E-02 

7.898E-02 

50 

3. 165E-03 

8.372E-02 

8.213E-02 

51 

3 .205E-03 

8 . 693E-02 

8.531E-02 

52 

3.240E-03 

9.017E-02 

8.854E-02 

53 

3.272E-03 

9 . 344E-02 

9 . 1 79E-02 

54 

3.31 3E-03 

9.675E-02 

9.508E-02 

55 

3.357E-03 

1 .0Q1E-01 

9.842E-02 

56 

3.390E-03 

1 .035E-01 

1 . 0 1 8E-0 1 

57 

3.400E-03 

1 .069E-01 

1 .052E-01 

58 

3.425E-03 

1 . 103E-01 

1 .086E-01 

59 

3.475E-03 

1 . 138E-01 

1 . 120E-01 

60 

3.500E-03 

1 . 173E-01 

1 . 155E-01 

61 

3.500E-03 

1 .208E-01 

1 . 1 90E“0 1 

62 

3.525E-03 

1 .243E-01 

1 .225E-0 1 

63 

3.575E-03 

1 .279E-01 

1 . 261E-0 1 

64 

3.625E-03 

1 .315E-01 

1 .297E-01 

65 

3.650E-0: 

1 .352E-01 

1 .333E-01 

66 

3.650E-03 

1 .388E-01 

1 . 370E-0 1 

67 

3 . 700E-03 

1 .425E-01 

1 .406E-0T 

68 

4.675E-03 

1 .472E-01 

1 444E-01 

69 

2.800E-03 

1 .500E-01 

1 . 500E-0 1 
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are about the same. The predicted results of mean axial velocity and turbulent kinetic energy at an axial 
location of 15 mm are shown in Figure 6.4.1-7. The results show that the turbulent kinetic energy has de- 
creased due to an excessive inlet dissipation rate in the inner region. On the other hand, the mean veloc- 
ity is not affected significantly. This can be attributed to the fact that the maximum value of the inlet e in 
the annular region has not been changed considerably. 

6.4.1 .2 The Algebraic Secon d-Moment Closure 

In this section, the predictions using the ASM have been compared with those from the k-e model. Simi- 
lar to the trends observed in the k-e model calculations, the effect of various discretization schemes on the 
ASM results was found to be rather insignificant. Consequently, results from different schemes will be 
shown only for some cases. The predicted mean axial velocity profiles from the ASM and k-e models 
have been compared in Figure 6.4.1-8. These results were obtained using the flux-spline discretization 
scheme. The use of ASM improves the overall agreement between the predictions and the experimental 
data. The major differences between the two turbulence models are seen in the regions where a maxima 
or a minima occurs in the velocity profiles. 

The predicted turbulent shear stress from the ASM has been compared with the experimental data in Fig- 
ure 6.4.1 -9. Here, results from the power-law differencing scheme have also been included to assess the 
numerical accuracy of the results. Both discretization schemes give nearly identical results, indicating 
that the solution is grid independent. The positive peak in the shear stress profile corresponds to the 
shear layer between the two streams and the negative peak corresponds to the shear layer associated with 
the expansion. The agreement between the calculation and the experimental data is good, although the 
peak values are not well predicted. 

The normal stresses at different axial locations are shown in Figure 6.4.1-10. Again, results from both the 
power-law and flux -spline schemes have been included. The two schemes perform essentially to the 
same level except for minor differences in the peak values. The differences between the calculated and 
experimental results are most significant in the case of the normal stresses. The experimental data indi- 
cate a faster mixing within the shear layer is predicted by the model. 

Comparisons of the radial and tangential Reynolds stress components with the measurements are pre- 
sented in Figures 6.4.1-11 and 6.4.1-12. Essentially, the same trends are observed for v 2 and w 2 profiles. 
The peak values of the w 2 , however, are reported higher. Such behavior has been correctly predicted by 
the ASM calculation. In addition, the maximum values are overpredicted for u 2 , underpredicted for v 2 , 
and closely predicted for w^. This clearly indicates the lack of performance of the pressure-strain model. 
One reason for this could be that the constant C 2 used is not suitable for complex turbulent flows. Since 
C 2 is determined from simple turbulent flows in local equilibrium, it would be more appropriate for 
equilibrium ASM than for nonequilibrium ASM. Another reason could be the incorrect modeling of the 
mean strain part of the pressure-strain term. It is very important to model the mean strain tensor 
correctly. Either the mean strain effects are not modeled at all or they should be accounted for properly 
(Nikjooy and So, 1987). Near the centerline region, substantial differences are observed between the 
model predictions and the experimental data. The discrepancies are attributed to inaccurate anisotropic 
diffusivity in momentum as well as k and e equations, an improper model used for the pressure-strain 
term, and the shortcoming of the algebraic stress model in predicting the normal stress components. 
Overall, the algebraic stress model can do a better job of predicting the shear stress components than the 
normal stresses. 

Numerical calculations were also performed starting from 15 mm downstream of the nozzle exit where 
no reverse flow occurred. This allowed the application of parabolic calculations with very fine grid to 
achieve numerically accurate solutions. The experimental profiles of mean and turbulence quantities at x 
= 15 mm were used as initial conditions for the parabolic calculations. The inlet turbulence dissipation 
rate profile was derived from Equation 69. 
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The k-e turbulence model was applied first for a single-phase coaxial jet calculation to establish the suit- 
ability of the model before introducing two-phase effects. The predicted mean axial velocity, turbulent 
kinetic energy, and shear stress profiles at the various axial stations, 15 (inlet plane), 25, 35, 50, 75, 150, 
and 300 mm, are shown in Figure 6.4.1-13 an< I have been normalized by the inner jet pipe radius and local 
centerline velocity. 

The k-e model yields good overall agreement with velocity measurements, however, the discrepancies are 
high for the kinetic energy distribution. The underpredicted values of the kinetic energy of turbulence 
are more obvious in the regions where the flow changes rapidly, in the developing shear layer between 
the two jets (first region) and where the shear layer reaches the centerline (second region). 


In the first region, the mixing layer between the two flows produced by the inner and outer streams is 
changing from a thin layer characterized by a high shear stress to a self-similar shear flow where the 
shear stress is relatively small. In both regions, the flow possesses two significant rate-of-strain compo- 
nents resulting from the velocity gradients in both the axial and radial directions. This type of flow is 
classified as complex shear flow, which is difficult to predict by using the standard version of the k-e 
model. 

The distribution of the turbulent shear stress and its dependence on the mean velocity gradient are also 
shown in Figure 6.4.1-13. According to the eddy-viscosity hypotheses, the sign of the shear stress and the 
velocity gradient is clearly consistent across the flow field. The positive shear stress in the inner jet region 
indicates that the flow still bears the characteristics of the upstream pipe flow, the gradient of mean axial 
velocity is still negative but approaching zero. The distribution of k and uv reach their more familiar 
shapes (as in a single free jet) when the mixing takes place across the entire flow and the mean axial ve- 
locity gradient becomes negative across the entire jet. 

The gas-phase rms velocities are shown in Figure 6.4.1-14. It is assumed that turbulent kinetic energy re- 
distributes with u 2 = k, v 2 = w 2 . The shapes of the three components have maximum values just behind 
the nozzle wall and minimum values at the central axis. But in the farther downstream region where the 
mixing layer merges with the central axis, the rms profiles take a uniform shape like that observed in fully 
developed jet flows. Figure 6.4.1-14 also shows that the radial rms velocity component, v 2 , has almost the 
same value of the axial component. 

The k-e model is extended to particle-laden jets' situations. The mass loading, defined as the ratio of par- 
tide-to-gas mass flow rate of the inner jet at the inlet plane, LR = 0.2, was considered. To distinguish be- 
tween the effects of mean and fluctuating gas velocity on particle transport, predictions were compared 
using ST and DT treatments to the measured data. 

Figure 6.4.1-15 presents the mean axial velocity, kinetic energy, and shear stress of the carrier phase. The 
high velocity of the external stream, compared to that of the central jet, causes a rapid increase in the axial 
velocity of the inner jet downstream of the exit plant. Because of this transfer of mass and momentum 
from tlW external to the internal stream, the carrier-phase velocity distribution where the particles exist 
becomes different from that of the axisymmetric jet flow case. This change in air velocity affects particle 
velocity, which is discussed in connection with Figure 6.4.1-15. 

The present particle loading is fairly small in terms of altering the carrier-phase flow field (Mostafa and 
Mongia, 1988). However, comparison of the last two planes in Figure 6.4.1-15 with those of Figure 6.4.1- 
13 shows some reductions in turbulence kinetic energy and shear stress in the near centerline region 
where the particles exist. This turbulence modulation is caused mainly by the fluctuating relative velocity 
between the particles and the carrier phase. Particles generally cause a reduction in gas turbulence and 
an increase in the dissipation rate of that energy. This phenomenon was simulated in the present study 
by introducing extra terms into the turbulence kinetic energy and its dissipation rate equations. Increas- 
ing the particle loading increases the turbulence modulation, as recently measured by Fleckhaus et al 
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(1987). Figure 6.4.1-15 shows that the present mathematical model yields fairly good agreement with the 
experimental data. 

Figure 6.4.1-16 presents the measured values of the three components of fluctuating velocity of the carrier 
phase when glass beads are present. The effect of the particles on the radial and azimuthal components is 
higher than that on the axial component. This condition may be attributed to the differences in the mean 
velocity strain profiles and to the small ratio of particle-to-gas velocity fluctuation. According to Fleck- 
haus et al (1987), the smaller the ratio of particle-to-gas velocity fluctuation, the higher is turbulence at- 
tenuation. This observation supports the postulate that particles increase the level of anisotropy and 
shows the need for detailed measurements of mean and fluctuating components of the two phases under 
different flow conditions. 

Figure 6.4.1-17 presents the data of the mean axial particle velocity. It also compares the predictions of 
the ST and DT treatments to the experimental data. The mean particle velocity profile is uniform over the 
entire cross section of the flow field. This behavior is different from that observed in round jet flow mea- 
surements and could be attributed to rapid mixing between inner and outer jets, which subject the parti- 
cles to a more uniform gas velocity distribution. This jet mixing process creates a strong negative radial 
velocity in the jet's outer region and causes the particle-number density to become narrower than the cor- 
responding profile for axisymmetric round jet flows. Figure 6.4.1-17 also shows that the ST provides 
good predictions compared to the experimental data, whereas the DT performs quite poorly at the 
downstream stations. That is, a particle moves radially because of its initial mean radial velocity and/ or 
the mean radial gas velocity, both of which are very small compared to the axial component. This effect 
explains the narrow distribution of particle mean axial velocity and number density predicted by DT. 

Figure 6.4.1-18 presents the measurements of the three components of fluctuating particle velocity and 
indicates the extent to which ST allows reasonable calculations. Although DT ignores these components 
entirely, ST simulates these particle-velocity fluctuations as a response to the carrier-phase components. 
The agreement between predictions and data is less than satisfactoiy at the last station, which could be 
attributed to the limitations of the assumptions embedded in the ST formation. Figure 6.4.1-18 also shows 
the high anisotropy of particle turbulent quantities that increase the anisotropy level of the carrier phase 
discussed in the analysis of Figure 6.4.1-16. 
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Figure 6.4. 1-1. Unconfined coaxial jets geometry. 
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Figure 6.4.1-2. Measured profiles at the inlet plane (x*2 mm). 
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Figure 6.4.1 -3. Unconfined coaxial jets grid layout. 
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Figure 6.4.1-4. Comparison of measurements with mean axial velocity calculations 

for unconfined coaxial jets. 


559 



o 



k (m 2 /s 2 ) 




k (nn 2 /s 2 } 


TE92-1251 


Figure 6.4.1-5. Comparison of measurements with turbulent kinetic energy calculations. 
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Figure 6.4.1 -6. Inlet dissipation rate profiles. 
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Figure 6.4. 1-7. Profiles of mean axial velocity and turbulent kinetic energy using different 

inlet turbulent dissipation rate. 
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Figure 6.4.1-8. Comparison of the k-c and the ASM predictions of axial velocity using 

flux-spline with measurements. 
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Figure 6.4.1-9. Comparison of the ASM prediction of turbulent shear stress with data. 
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Figure 6.4.1-10. Comparison of the ASM prediction of axial normal stress with data. 
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Figure 6.4.1-12. Comparison of the ASM prediction of tangential normal stress with data 
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Figure 6.4.1-13. Radial profiles of normalized mean axial velocity, turbulent kinetic energy, and shear 

stress for the single-phase coaxial jets. 
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Figure 6.4.1-15. Radial profiles of normalized mean axial velocity, turbulent kinetic energy, and shear 

stress for the particle-laden coaxial jets. 
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Figure 6.4.1-17. Radial profiles of normalized particle axial velocity. 
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In this section comparison of computational results with data is presented for coaxial flow with and with- 
out particle beads. For the present study, a confined flow configuration was selected. In this configura- 
tion, the injector was directed vertically downward. The experimental test parameters are summarized in 
Table 6.4-III, and the sketch of the test section is shown in Figure 6.4.2-1. 

The set of governing partial differential equations apply for nonswirling continuous (carrier) phase flow 
consist of equations for continuity, axial (x) and radial (r) momenta, e, and four Reynolds stress compo- 
nents. The finite volume approach (Patankar, 1980) is used to reduce the continuous equations to a set of 
coupled discrete equations. The numerical solutions are obtained using the flux-spline differencing 
scheme (FSDS) (Varejao, 1979). 

The elliptic nature of transfer equations requires that boundary conditions be specified on the four sides 
of the solution domain. Tourlclnds of boundaries, inlet, axis of symmetry, wall, and the outlet, need con- 
sideration. At the inlet boundary, which is located at the first measurement plane (x = 4 mm), the mea- 
sured profiles of U, V, W, u^, v^, w^, and uv are applied. The inlet dissipation rate is prescribed based on 
the assumption of constant length scale and the turbulent kinetic energy, namely 


0.2 Ra 

where R3 (R3 = 76 mm) is the pipe radius. The macrolength scale of 0.2 R3 was estimated through the 
sensitivity analysis. 

At the axis of symmetry, the radial velocity, shear stresses, and the radial gradients of other variables are 
set to zero. All stream wise gradients are presumed zero in the exit plane of the calculation domain except 
for axial velocity which needs to be corrected to satisfy integral mass balance. The conventional loga- 
rithmic law, which is based on the local equilibrium assumption, is applied to determine the wall shear 
stress component The shear stress is then used as the boundary condition for U momentum and uv 
equations. The boundary conditions for normal stresses are imposed through the modifications of their 
production terms using the new calculated wall shear stress value. 


Table 6.4-m. 

nfined coaxial jets experimental conditions at 4 mm downstream of nozzle exit. 


Continuous phase 

Medium 

Air 

Density, r 

1.17 Kg/m 3 

Inner jet mass flow rate, m] 

0.0021 Kg/s 

Annular jet mass flow rate, m2 

0.0033 kg/s 

Dilute jet mass flow rate, m3 

0.0272 kg/s 

Averaged velocity of inner jet, U Q 

3.935 m/ s 

Discrete phase 

Medium 

Glass beads 

Density, rp 

2500 kg/m 3 

Beads diameter. Dp 

105 mm 

Centerline velocity 

4.2 m/s 

Centerline beads rate 

901/s 

Mass flow rate 

0.00193 kg/s 

Mass loading beads to air, LR 

0.925 
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The coupled equations and boundary conditions are solved numerically in a sequential manner using the 
staggered grids for velocities and shear stresses. The main advantage of staggering the locations of 
stresses is to enhance the numerical stability, a result of high coupling between the shear stresses and re- 
lated mean strains. The iteration sequence employs the SIMPLER algorithm (Patankar, 1980) to handle 
the coupling between the continuity and momentum equations. The algebraic equations are solved using 
a line-by-line TDMA. 

The ordinary differential equations describing the dispersed phase are solved using forward numerical 
integration. The calculations start from the first measurement plane (x = 4 mm) where the mean and 
fluctuation velocity profiles of particles and their rates of injection are available. At the wall, particles are 
assumed to bounce back with an angle of reflection equal to the angle of incidence. 

Numerical iterations are performed over the continuous and dispersed phases until the absolute sums of 
the normalized mass and momentum of the carrier phase at all internal grid points, and the change of the 

particle source terms, are less than 10'^. 

A nonuniform grid of 61 x 57 in the axial and radial directions is used for continuous phase calculations 
(Figure 6.4.2-2). A finer grid sparing is assigned near the inlet, centerline, and in the shear layer. The 
computational domain extends from the first experimental location (x = 4 mm) to 300 mm downstream of 
the jet exit. Since the measured flow does not show any x dependence at x > 250 mm, the specified condi- 
tion 3 /3x = 0 at the exit plane of the calculation domain is reasonable. The discrete phase calculation also 
starts from the same axial location, i.e., x = 4 mm. Four thousand particles distributed at forty radial loca- 
tions are used for ST treatment. 

The present calculations have been made using the flux-spline differencing scheme and the calculated re- 
sults are essentially free of numerical diffusion. Therefore, the discrepancy between the experimental 
data and the prediction can be attributed to two sources, improper boundary conditions at the inlet plane 
and the deficiency of the turbulence model. As regards the inlet conditions, all quantities except the dis- 
sipation rate (e) were prescribed from the measurement. These profiles are shown in Figure 6.4.2-3 and 
have been normalized by the inner jet averaged velocity (Uq) and the pipe radius (R3)- Above r/R3>0.5, 
the flow field is similar to plug flow condition created by the strong coflow. This region is not the focus 
of the study and has not been shown in the presented figures. The uncertainties in the derivation of inlet 
e profile would adversely affect the calculation at downstream locations. It was shown, however, that the 
use of inlet e derived from the constant length-scale assumption can result in a better prediction (Nikjooy 
and Mongia, 1991). 

Comparison of the normalized axial mean velocity profiles for single-phase flow with experimental data 
along the streamwise direction is presented in Figure 6.4.2-4. These results were obtained from the k-e 
model. The predictions are in very good agreement with measured values. There are some discrepancies 
near the maximum flow region which can be due to the higher radial diffusion process. The flow dis- 
charging from the inner pipe has a flat velocity profile which decelerates gradually due to the pressure ef- 
fects. Comparisons of the turbulent shear stress and kinetic energy are shown in Figures 6.4.2-S and 6.4.2- 
6. The k-e model overestimated the maximum level of the turbulent kinetic energy near the nozzle exit, 
however, the agreement between the calculated shear stress and data is excellent. It seems that the over- 
predicting of turbulent kinetic energy in the upstream region resulted in higher effective diffusion, which 
is responsible for underestimating the peak velocities. At about x = 150 mm both flow streams have been 
completely mixed. The agreement between the predicted and experimental values of kinetic energy is not 
as good as that for the axial velocity. Even though the trends are similar, the predicted kinetic energy is 
smaller than that derived from the measurements in downstream regions. 

Comparison of the mean and turbulence quantities by the DSM closure with the measured values is pre- 
sented in Figures 6.4.2-7 through 6.4.2-11. The predicted mean velocity profiles are in excellent agreement 
with the data. The peak and centerline velocities were very well predicted. The pressure-strain correla- 
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tion model used for these calculations is just based on the retum-to-isotropy model (Cl = 5.0, C 2 = 0.0). 
The differences between the calculated and experimental results are most significant in the case of the 
normal stresses. The maximum axial turbulence intensity has been overpredicted in most of the flow re- 
gion, however, the trend was accurately predicted. This behavior resulted in underprediction of the 
maximum values of other rms components, especially the radial fluctuation. Essentially, the same trends 
are observed for radial and tangential stress components. The discrepancies between the model and data 
clearly indicate the lack of performance of the pressure-strain model. One reason for this could be the 
lack of an appropriate model for the "rapid part" which was not considered here. The turbulent shear 
stress was accurately predicted. The negative peak in the shear stress profile corresponds to the shear 
layer between the two streams and the positive peak corresponds to the shear layer associated with the 
expansion. 

Comparisons of the normalized mean axial velocity and Reynolds stress predicted by DSM and ASM clo- 
sures with the measured values are presented in Figures 6.4.2-12 through 6.4.2-16. The calculated velocity 
fields predicted by both models are very similar at downstream locations. Some differences are observed 
between the two models in the maximum an J expansion zones at upstream regions. Overall, the ASM 
prediction is in good agreement with data. Examination of the calculated rms profiles indicates that the 
DSM’s results mimic the experimental data better at all stations. A wide disparity exists between the 
measured normal Reynolds stresses and the ASM's results, especially in the peak zone. As regards the 
turbulent shear stress, the ASM underpredicted the maximum and minimum data points in upstream lo- 
cations. The predicted turbulent kinetic energy by the k-e, DSM, and ASM (Figure 6.4.2-17) shows that 
the DSM simulated the data letter than other models. 

The calculation for a particle-laden jet was also performed over the same computational domain and the 
same grid distribution was used. The inlet conditions, except the dissipation rate (e) were specified from 
the experimental data. These profiles are shown in Figure 6.4.2-18 and have been normalized by the inner 
jet averaged velocity (U 0 ), particles flux at the centerline (N 0 ), and the pipe radius (R 3 ). The inlet dissipa- 
tion rate is prescribed based on the assumption of constant length scale and the turbulent kinetic energy. 

The normalized axial component of mean velocity by DSM at different streamwise locations is shown in 
Figure 6.4.2-19. The predicted result is in excellent agreement with die experimental data. The velocity 
profiles are quite flat near the centerline. The measurement shows that the centerline velocity initially de- 
creases about 10% from the inlet plan to station x = 25 mm. The deceleration in axial velocity probably 
results from the pressure effect. The predicted pressure distribution shows a negative radial gradient fa- 
voring an outward motion of the flow. As a result, the axial velocity is slowing down to satisfy the mass 
conservation. The centerline velocity at subsequent streamwise locations then increases, indicating the 
momentum transfer from the annular jet to the core region. It should be mentioned that the choice of 
model constants in the pressure redistribution term has been found to be effective especially on the turbu- 
lence field (Nikjooy and Mongia, 1991). The case of coaxial jets was also calculated for Cl = 1.5 and C 2 = 
0.6 in the <fcj model. The prediction however was in large discrepancy with the data. 

Comparisons of the normalized rms velocity components predicted by DSM closure with the measured 
values are presented in Figures 6.4.2-20, 6.4.2-21, and 6.4.2-22. Near the injector exit, e.g. x<15 mm, the 
axial rms peaks at two radial locations correspond to the sharp decrease in the mean velocity. The pro- 
duction of turbulent kinetic energy is increased at these locations, consistent with the maximum turbulent 
intensity profiles observed in the measurement. Overall, the calculated Vu^ profiles are in good agree- 
ment with the data. The peak values are slightly overpredicted at some locations. This behavior resulted 
in underprediction of the maximum values of other rms components, especially the radial fluctuation. 

The discrepancy observed in the turbulence field between the prediction and the measurement is related 
in large to the deficiency of the pressure-strain model. The turbulence quantities however were all well 
predicted near the centerline. At downstream locations, turbulent intensities are relatively uniform 
across the pipe, indicating a well-mixed flow condition. As regards the turbulent shear stress, uv, the cal- 
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culated profiles are very similar to the data behavior (Figure 6.4.2-23). The results demonstrate that the 
dispersed phase had no significant effects on the continuous phase. 

Predicted and measured mean particle velocities along the axis are shown in Figures 6.42-24 and 6.4.2-25. 
The calculated velocities are based on the particles' mass flow rate weighted averages. Only a portion of 
the entire radial section where particles could be found are presented. The predicted axial values are in 
good agreement with data. Despite the variations observed in the gas phase velocity near the centerline, 
the particle velocity shows almost no change. This is related to the fact that the particle dynamic relax- 
ation time (x<j) is very large compared to the turbulent characteristic time (x). The reason for large x<i is 
found in the large particle-to-gas density ratio (2500:1.17). This clearly shows that the inertia force is re- 
sponsible for the particles' motion and the effects of drag force are marginal. The predicted profiles of 
mean radial velocity of particles are in qualitative agreement with the data. The radial velocity compo- 
nent however is at least one order of magnitude smaller than the axial component. It is therefore reason- 
able to believe that their quantitative disagreement would not severely affect other results. 

The agreements between the predicted and measured fluctuating velocities of the particles are not favor- 
able in downstream locations (Figures 6.4.2-26, 6.4.2-27, and 6.4.2-28). The data show similar radial and 
tangential fluctuation components. On the other hand, the axial rms value is about 15 times larger than 
others. The model has successfully predicted the anisotropy feature of the dispersed flow field. How- 
ever, all three components are predicted to be lower than the data. The results for the normalized particle 
number density are illustrated in Figure 6.4.2-29. The reference value used for normalization is the parti- 
cle number density on the centerline at the inlet plane (x = 4 mm). The predicted profiles presented here 
also represent the particles' mass flow rate per unit area or mass flux. The predicted mass flux is in rea- 
sonable agreement with data at upstream locations (x < 35 mm). However, the calculation shows a larger 
spread of the particles at the downstream region. As a result, the particles mass flux is underpredicted in 
the central region of the flow at those places. Although the gas phase was well-mixed downstream of the 
flow field, the particles will tend to accumulate near the centerline. The reasons can be due to the lack of 
lift force, which was discarded from the interfacial force term. 

In the present study, the Reynolds equations have been extended to calculate two-phase turbulent flow 
without evaporation. The selected model used in this study is the extension of those desired for constant 
density single-phase turbulent flow. Therefore, some of the discrepancies observed in the results are re- 
lated to this matter. However, the errors arising from the inlet dissipation rate profile and/ or the model 
constants used in the pressure redistribution term should not be overlooked. Overall, the mean and tur- 
bulence fields associated with the continuous phase were reasonably well predicted by DSM. The gas 
appeared to be unaffected by the particles' motion. The discrete phase results are also encouraging. The 
discrepancies observed in the fluctuation velocities and the mass flux profiles for the particles are more 
clear in the downstream region where the carrier phase flow is well mixed. However, the particles still 
tend to accumulate near the centerline region. This shortcoming may be removed by including the effects 
of the lift force into the carrier phase equations. 
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Figure 6.42-2. Confined coaxial jets - grid layout. 
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Figure 6.4.2-5. Comparison of predicted turbulent shear stress by k-e model with data. 










Figure 6.42-6. Comparison of predicted turbulent kinetic energy by k-c model with data. 



Figure 6.4.2-7. Comparison of predicted mean axial velocity by DSM with data. 




Figure 6.4.2-8. Comparison of predicted axial rms velocity by DSM with data. 





.4.2-9. Comparison of predicted radial rms velocity by DSM with data. 




Figure 6.4.2-10. Comparison of predicted tangential rms velocity by DSM with data. 











Figure 6.42-12. Comparison of predicted mean axial velocity by DSM and ASM with data. 





Figure 6.4.2-14. Comparison of predicted radial rms velocity by DSM and ASM with data. 
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Figure 6.42-15. Comparison of predicted tangential rms velocity by DSM and ASM with data. 




Figure 6.4.2-16. Comparison of predicted turbulent shear stress profiles by DSM and ASM with data. 




figure 6.4.2-17. Comparison of predicted turbulent kinetic energy by k-e, DSM, and ASM with date 






continuous phase by DSM with data. 



Figure 6.4.2-20. Comparison of calculated normalized stream wise turbulence intensity of 

continuous phase by DSM with data. 
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dispersed phase by DSM with data. 
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Figure 6.4.2-28. Comparison of calculated normalized tangential turbulence intensity of 

dispersed phase by DSM with data. 




6.5 UNCONFINED COAXIAL JETS WITH SWIRLING ANNULAR FLOW 


In this section, computations for unconfined axisymmetric swirling coaxial jets are reported. In this con- 
figuration, an axial jet injector surrounded by a swirling annular jet, as shown in Figure 6.5-1, was di- 
rected vertically downward within a 457 mm^ wire mesh screen. The entire test assembly was sur- 
rounded by a flexible plastic enclosure. This helped damp-out extraneous room drafts, and more impor- 
tantly, allowed uniform seeding of the entrained air, thereby permitting unbiased measurements in the 
outer region. Data were obtained at eight axial stations 4, 15, 25, 35, 50, 75, 100, and 150 mm from the exit 
plane of the injector. At each axial station, between 10 and 20 radial points were scanned. The flow con- 
ditions used for this case are given in Table 6.5-1. The center jet (diameter, D = 24.1 mm) is surrounded by 
an annular jet with the inner and outer diameters of 29 and 36.7 mm, respectively. The effective area and 
axial velocity ratios of the annular jet to the center jet are 0.87 and 1.8, respectively. 

The control-volume approach was applied to discretize the governing equations. To minimize the effects 
of false diffusion, numerical solutions have been obtained using a flux-spline scheme for the convection 
and diffusion terms in the transport equations on a fine grid. The computational mesh used for all calcu- 
lations consisted of 76 x 69 nonuniformly distributed grid points in the axial and radial directions. A finer 
grid spacing was used near the inlet and the centerline in the shear layer and within the recirculation 
zone (Figure 6.5-2). The tabulated axial (x) and radial (r) grid points are shown in Table 65-11. The com- 
putational domain extended from the first measurement plane, located downstream of the nozzle exit at a 
distance of 3.0 mm, to 40 inner jet diameters downstream of the nozzle exit. 

A calculation procedure for elliptic flow requires boundary conditions on all boundaries of the computa- 
tional domain. Four kinds of boundaries, namely inlet, axis of symmetry, outlet, and the entrainment 
boundary, need consideration. At the inlet boundary, which was located at the first measurement plane, 
the measured profiles of mean velocity and Reynolds stress components were prescribed. These profiles 
are shown in Figure 65-3. The k-profile was obtained from the measured Reynolds stresses. This kinetic 
energy distribution was used to derive the e values at the inlet plane through the relationship 

e = K 15 (71) 

0.2 Dj 

where Dj is the inner jet diameter. This approach is based on the assumption of constant length scale. 

At the axis of symmetry, shear stresses and radial gradients of other variables are set to zero. At the out- 
let, axial diffusion is neglected for all variables. In the radial direction, the entrainment boundary was 
placed at a distance of six jet diameters from the axis of symmetry. Along the entrainment boundary, the 
quantity (rV) was assumed constant. In addition, the axial velocity U and tangential velocity W were as- 
sumed zero. Turbulent kinetic energy and its dissipation rate were assigned arbitrarily low values yield- 
ing an eddy viscosity, m = 10 p. 


Table 6.5-1. 

Unconfined swirling coaxial jets experimental conditions. 


Average inner velocity, m/s 

3.9 

Average annular velocity, m/ s 

7.1 

Centerline velocity U 0 m/ s 

3.93 

Density r, kg/m 3 

1.178 

Inner jet mass flow rate, mi, kg/m 3 

0.0021 

Annular jet mass flow rate, m2, kg/m 3 

0.0033 

Reynolds number. Re = 4ml / pmD 

5712 

Swirl angle, deg 

60 
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Tabic 6.5-11. 

Unconfined swirling coaxial jets grid d efinition^ 


STREAMWISE COORDINATES OF THE GRID 


I 

DX 

X 

XU 

1 

0 . 000E+00 

0 .OOOE+OO 

0. 000E+00 

2 

7 . 750E-04 

7.750E-04 

0. OOOE+OO 

3 

1 .625E-03 

2 - 400E-03 

1 .550E-03 

4 

2.050E-03 

4.450E-03 

3 . 250E-03 

5 

2.500E-03 

6 . 950E-03 

5 . 650E-03 

6 

2 . 700E-03 

9.650E-03 

8 , 250E-03 

7 

2.938E-03 

1 . 259E-02 

1 .105E-02 

8 

3.313E-03 

1 . 590E-02 

1 .413E-02 

9 

3.88BE-03 

1 . 979E-02 

1 . 768E-02 

10 

4.58BE-03 

2 . 438E-02 

2. 190E-02 

11 

5.250E-03 

2 . 963E-02 

2.685E-02 

12 

5.775E-03 

3.540E-02 

3.240E-02 

13 

6.23BE-03 

4. 164E-02 

3.840E-02 

14 

6 . 713E-03 

4 . 835E-02 

4 . 488E-02 

15 

7 . 063E-03 

5.541E-02 

5.183E-02 

16 

7.200E-03 

6.261E-02 

5.900E-02 

17 

7 . 250E-03 

6.986E-02 

6.623E-02 

18 

7 . 28BE-03 

7.715E-02 

7.350E-02 

19 

7.300E-03 

8 . 445E-02 

8.080E-02 

20 

7 . 300E-03 

9.175E-02 

8.810E-02 

21 

7.288E-03 

9 . 904E-02 

9 . 540E-02 

22 

7.275E-03 

1 .063E-01 

1 .027E-01 

23 

7.288E-03 

1 . 136E-01 

1 . 100E-01 

24 

7.300E-03 

1 . 209E-0 1 

1 . 173E-01 

25 

7.300E-03 

1 . 282E-0 1 

1 . 246E-01 

26 

7 . 300E-03 

1 , 355E-0 1 

1 .319E-01 

27 

7.300E-03 

1 . 428E-0 1 

1 .392E-01 

28 

7 . 300E-03 

1 .501E-01 

1 . 465E-0 1 

29 

7.300E-03 

1 .574E-01 

1 .538E-01 

30 

7.300E-03 

1 .647E-01 

1 .611E-01 

31 

7.300E-03 

1 . 720E-0 1 

1 .684E-01 

32 

7.300E-03 

1 .793E-01 

1 .757E-01 

33 

7.300E-03 

1 .866E-01 

1 . 830E-0 1 

34 

7.300E-03 

1 .939E-01 

1 . 903E-0 1 

35 

7.300E-03 

2.012E-01 

1 .976E-01 

36 

7.300E-03 

2.085E-01 

2.049E-01 

37 

7 . 300E-03 

2. 158E-01 

2.122E-01 

38 

7 . 300E-03 

2.231E-01 

2. 195E-01 

39 

7.200E-03 

2.303E-01 

2.268E-01 

40 

7.200E-03 

2.375E-01 

2.339E-01 

41 

7.400E-03 

2.449E-01 

2.412E-01 

42 

7.400E-03 

2.523E-01 

2.487E-01 

43 

7.325E-03 

2.596E-01 

2 . 560E-0 1 

44 

7 . 325E-03 

2.670E-01 

2.633E-01 

45 

7.275E-03 

2.742E-01 

2.706E-01 

46 

7.275E-03 

2.815E-01 

2 . 779E-0 1 

47 

7.300E-03 

2.888E-01 

2.852E-01 

48 

7 . 450E-03 

2.963E-01 

2 . 925E-0 1 
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Table 63-11 (cont). 


49 

7.775E-03 

3.040E-01 

3 . 001E-01 

50 

7.975E-03 

3. 120E-01 

3 . 080E-01 

51 

8.500E-03 

3.205E-01 

3 . 160E-01 

52 

9.500E-03 

3 . 300E-01 

3.250E-01 

53 

1 . 000E-02 

3 . 400E-01 

3.350E-01 

54 

1 . OOOE-02 

3 . 500E-0 1 

3 , 450E-0 1 

55 

1 . 125E-02 

3.613E-01 

3.550E-01 

56 

1 . 375E-02 

3 . 750E-01 

3.675E-01 

57 

1 .500E-02 

5. 9001-61 

3 . 825E-01 

58 

1 .500E-02 

4.050E-01 

3.975E-01 

59 

1 .500E-02 

4.200E-01 

4.125E-01 

60 

1 .500E-02 

4.350E-01 

4.275E-01 

61 

7.500E-03 

4.425E-01 

4.425E-01 


TRANSVERS 

COORDINATES OF 

THE GRID 


J 

DY 

Y 

YV 

1 

0 . 000E+00 

5.000E-04 

0 . 000E+00 

2 

2.500E-04 

7.500E-04 

5 . 000E-04 

3 

7.500E-04 

1 :500l>03 

1 . 000E-03 

4 

1 . 00GE-O3 

2.500E-03 

2. 000E-03 

5 

1 . OOOE-03 

3.500E-03 

3. 000E-03 

6 

1 . 000E-03 

4 ; 500E-03 

4. OOOE-03 

7 

1 .000E-03 

5.500E-03 

5. 000E-03 

8 

1 . 000E-03 

6.500E-03 

6. 000E-03 

9 

1 . 000E-03 

7.500E-03 

7. 000E-03 

10 

1 . 000E-03 

8 , 500E-03 

8. 000E-03 

1 1 

1 . 000E-03 

9.500E-03 

9. 000E-03 

12 

1 . 000E-03 

1 . 050E-02 

1 .000E-02 

13 

1 . 000E-03 

1 . 150E-02 

1 .100E-02 

14 

1 . OOOE-03 

1 . 250E-02 

1 . 200E-02 

15 

7.500E-04 

1 . 325E-02 

1 . 300E-02 

16 

5 . 000E-04 

1 , 375E-02 

1 . 350E-02 

17 

5 . 0Q0E-04 

1 . 425E-02 

1 . 400E-02 

18 

5.000E-04 

1 . 475E-02 

1 .450E-02 

19 

5.000E-04 

1 .525E-02 

1.500E-02 

20 

5.000E-04 

1 .575E-02 

1 . 550E-02 

21 

5.000E-04 

1 . 625E-02 

1 . 600E-02 

22 

5.000E-04 

1 .675E-02 

1 .650E-02 

23 

5.000E-04 

1 .725E-02 

1 . 700E-02 

24 

5.000E-04 

1 .775E-02 

1 . 750E-02 

25 

5.000E-04 

1 .825E-02 

1 . 800E-02 

26 

5.000E-04 

1 .875E-02 

1 . 850E-02 

27 

5.000E-04 

1 . 925E-02 

1 . 900E-02 

28 

5.000E-04 

1 .975E-02 

1 .950E-02 

29 

7.500E-04 

2.050E-02 

2. 000E-02 

30 

l'. 000E-03 

2. 150E-02 

2.100E-02 

31 

1 . OOOE-03 

2 . 250E-02 

2.200E-02 

32 

1 . 000E-03 

2.350E-02 

2 . 300E-02 


PART 2 OF 3 
TE92-12M 


610 



Table 6.5-II (cont). 


33 

1 .000E-03 

2.450E-02 

2.400E-02 

34 

1 . OOOE-03 

2.550E-02 

2.500E-02 

35 

1 . 000E-03 

2.650E-02 

2 . 600E-02 

36 

1 . 000E-03 

2 . 750E-02 

2 . 700E-02 

37 

1 .000E-03 

2.850E-02 

2.800E-02 

38 

1 .OOOE-03 

2.950E-02 

2 . 900E-02 

39 

1 .500E-03 

3. 100E-02 

3.000E-02 

40 

2. 000E-03 

3.300E-02 

3 . 200E-02 

41 

2.500E-03 

3.550E-Q2 

3.400E-02 

42 

3. 000E-03 

3.850E-02 

3.700E-02 

43 

3. 000E-03 

4.150E-02 

4.000E-02 

44 

3. 000E-03 

4.450E-02 

4 . 300E-02 

45 

3.500E-03 

4.800E-02 

4 . 600E-02 

46 

4. 000E-03 

5.200E-02 

5.000E-02 

47 

4.500E-03 

5.650E-02 

5.400E-02 

48 

5. 000E-03 

6.150E-02 

5.900E-02 

49 

5. 000E-03 

6.650E-02 

6.400E-02 

50 

5. 000E-03 

7.150E-02 

6.900E-02 

51 

5. OOOE-03 

7.650E-02 

7.400E-02 

52 

5. 000E-03 

8. 150E-02 

7.900E-02 

53 

5. 000E-03 

8 . 650E-02 

8.400E-02 

54 

5. 000E-03 

9. 150E-02 

8.900E-02 

55 

5. 000E-03 

9.650E-02 

9.400E-02 

56 

5.500E-03 

1 . 020E-01 

9.900E-02 

57 

5.500E-03 

1 .075E-01 

1 .050E-01 

58 

5. 000E-03 

1 . 125E-01 

1 . 100E-01 

59 

5. 000E-03 

1 . 175E-01 

1 .150E-01 

60 

5. OOOE-03 

1 .225E-01 

1 .200E-01 

61 

5. 000E-03 

1 .275E-01 

1 .250E-01 

62 

5. 000E-03 

1 .325E-01 

1 . 300E-01 

63 

5. 000E-03 

1 .375E-01 

1 .350E-01 

64 

5. 000E-03 

1 .425E-01 

1 .400E-01 

65 

5 . 000E-03 

1 .475E-01 

1 .450E-01 

66 

2.500E-03 

1 .500E-01 

1 .50GE-01 
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Predictions of mean and turbulence fields obtained from the k-e model are discussed first. Comparison of 
the calculated mean axial and tangential (azimuthal) velocity profiles with the experimental data is pre- 
sented in Figures 6.5-4 and 6.5-5. In this case, the k-e model fails to display the size and strength of the 
experimental recirculation zone. The k-e model mimics the data trend reasonably well away from the 
symmetry axis, however, the peaks were slightly underpredicted. The underprediction for a centerline 
velocity stems from the incorrect representation of the turbulent diffusion process. The radial normal 
stress is particularly important in the upstream region and, as a consequence, the isotropic viscosity hy- 
pothesis is inadequate. The k-e model prediction of tangential velocity is in relatively good agreement 
with data. The calculated profiles are in accordance with the data trend, although the maximum values 
are slightly overpredicted. At the farthest downstream locations, the measured tangential profiles are 
skewed toward the centerline, however, the k-e model fails to demonstrate such a behavior. 

Comparison of the predicted turbulent kinetic energy, k, with the data is shown in Figure 6.5-6. The 
agreement between the predicted and experimental values of k is not as good as that for the velocity. The 
discrepancies near the centerline region can be attributed to the lack of well-predicted velocity field, inad- 
equacy of the eddy viscosity hypothesis, and turbulent kinetic dissipation rate, e. The inlet distribution of 
length scale has a significant effect on the e distribution in the upstream region where the diffusion terms 
have a large effect on the mean velocity, it has no significant effect in the downstream region where the 
diffusion effects are small. The k-e model tends to connect the dissipation rate too strongly to the local 
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mean velocity field (Habib and Whitelaw, 1979) which is inappropriate since the dissipation occurs in the 
finest scales of motion, and these do not reflect the local mean strain field. The predicted turbulent shear 
stress, uv, profiles have been compared with the experimental data in Figure 65-7. The negative peak in 
the shear stress profile corresponds to the shear layer between the two streams and the positive peak cor- 
responds to the shear layer associated with the expansion. The agreement between the calculation and 
the experimental data is good, although the peak values are not well predicted. 

Having demonstrated the results of the k-e model, attention is now turned to the performance of the DSM 
closure. Results of the mean velocity components obtained from the k-e model and the DSM are com- 
pared with the experimental data (Figures 65-8 and 65-9). A wide disparity exists between the models' 
prediction of the axial velocity near the centerline. It is evident that there are substantial differences in 
the capability of the various models to promote or hinder formation of the recirculation zone. It is clear 
that the proper turbulence model is dependent on the location within the flow field. With regard to the 
comparison between measurements and calculations, the predictions by DSM seem to be slightly better 
for tangential velocity and the central recirculation zone, whereas those by k-e model are closer to exper- 
imental data for maximum velocity. An examination of the calculated Reynolds stresses (Figures 65-10, 
6.5-11, and 6.5-12) indicates that the relative performance of the model is dependent on the flow region. 
The results show differences especially in the prediction of the centerline values in the developing region. 

The reason can be attributed to the inlet dissipation rate distribution. Numerical study indicates that the 
inlet e profile is a very important factor in predicting the maximum values of mean and turbulence quan- 
tities. In general, the DSM closure maintains the correct shape of Reynolds stress profiles which are con- 
sistent with the behavior of the mean velocity field. 

Comparisons of calculated shear stresses with data are shown in Figures 6.5-13 and 6.5-14 . In compari- 
son with the data, the predicted values are in good agreement in the developed region, however, the loca- 
tions of the maxima and minima have been shifted in the developing region of the flow field. In the outer 
part of the recirculation zone, where the effect of entrainment is more pronounced, and near the axis of 
symmetry, where the velocity gradient is negligible in the radial direction, the discrepancies between the 
k-e model and DSM results are in fact diminished. 
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Figure 6.5-1. Unconfined axisymmetric swirling coaxial jets geometry. 


613 



150 mm 



>s 


II 

>N 


614 


Figure 6-5-2. Unconfined swirling coaxial jets grid layout 



Figure 65-3. Measured profiles at the inlet plane (x = 4 mm) (1 of 2). 








Figure 65-3. Measured profiles at the inlet plane (x = 4 mm) (2 of 2). 
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Figure 65-4. Comparison of calculated mean axial velocity by k-e model with data (1 of 4). 
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Figure 65-4. Comparison of calculated mean axial velocity by k-e model with data (2 of 4). 
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Figure 65-5. Comparison of calculated mean tangential velocity by k-e model with data (1 of 4). 
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Figure 65-5. Comparison of calculated mean tangential velocity by k-e model with data (2 of 4). 









Figure 65-6. Comparison of predicted turbulent kinetic energy by k-e model with measurement (1 of 4). 



Figure 6.5-6. Comparison of predicted turbulent kinetic energy by k-e model with measurement (2 of 4). 


RvmmJ R(mm) 

30.0 60.0 90.0 <§' 0.0 20.0 40.0 



x = 50. mm 




x = 75. mm 


K(m2/s2) 


2.0 

K(m2/s2) 


TE92-1311 


ire 65-6. Comparison of predicted turbulent kinetic energy by k-e model with measurement (3 of 4). 
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Figure 65-6. Comparison of predicted turbulent kinetic energy by k-e model with measurement (4 of 4). 
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Figure 6.5-7. Comparison of predicted turbulent shear stress by k-e model with data (1 of 4). 
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Figure 6.5-7. Comparison of predicted turbulent sh<ar stress by k-e model with data (2 of 4). 
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Figure 6.5-7. Comparison of predicted turbulent shear stress by k-e model with data (3 of 4). 
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Figure 6.5-7. Comparison of predicted turbulent shear stress by k-e model with data (4 of 4). 
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Figure 65-8. Comparison of calculated mean axial velocity by k-e and DSM with data (1 of 4). 
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Figure 65-8. Comparison of calculated mean axial velocity by k-e and DSM with data (2 of 4). 
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Figure 65-9. Comparison of calculated mean tangential velocity by k-e and DSM with data (1 of 4 ). 



Figure 65-9. Comparison of calculated mean tangential velocity by k-e and DSM with data (2 of 4). 
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jure 65-9. Comparison of calculated mean tangential velocity by k-e and DSM with data (3 of 4). 
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Figure 65-9. Comparison of calculated mean tangential velocity by k-e and DSM with data (4 of 4). 
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are 6.5-10. Comparison of calculated axial rms velocity component by DSM with data (1 of 4). 



Figure 6.5-10. Comparison of calculated axial rms velocity component by DSM with data (2 of 4). 





Figure 6.5-10. Comparison of calculated axial rms velocity component by DSM with data (3 of 4). 
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Figure 6.5-10. Comparison of calculated axial rms velocity component by DSM with data (4 of 4). 
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jure 65-11. Comparison of calculated radial rms velocity component by DSM with data (1 of 4) 



Figure 65-11. Comparison of calculated radial rms velocity component by DSM with data (2 of 4). 
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Figure 63-1 1 . Comparison of calculated radial rms velocity component by DSM with data (3 of 4). 



Figure 63-1 1 . Comparison of calculated radial rms velocity component by DSM with data (4 of 4). 
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Figure 6.5-12. Comparison of calculated tangential rms velocity component by DSM with data (1 of 4). 



Figure 6.5-12. Comparison of calculated tangential rms velocity component by DSM with data (2 of 4). 
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Figure 6.5-12. Comparison of calculated tangential rms velocity component by DSM with data (3 of 4). 



Figure 6.5-12. Comparison of calculated tangential rms velocity component by DSM with data (4 of 4). 
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Figure 6.5-13. Comparison of calculated uv profiles by k-e and DSM with data (1 of 4). 
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Figure 6.5-13. Comparison of calculated uv profiles by k-e and DSM with data (2 of 4). 
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Figure 6.5-13. Comparison of calculated uv profiles by k-e and DSM with data (3 of 4). 
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Figure 6.5-13. Comparison of calculated uv profiles by k-c and DSM with data (4 of 4). 
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Figure 65-14. Comparison of calculated uw profiles by DSM with data (1 of 4). 



Figure 65-14. Comparison of calculated uw profiles by DSM with data (2 of 4). 




Figure 6.5-14. Comparison of calculated uw profiles by D5M with data (3 of 4). 
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Figure 65-14. Comparison of calculated uw profiles by DSM with data (4 of 4). 
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6.6 AIRBLAST INJECTOR 


In this section, computations for the single-phase flow exhausted from the airblast injector are reported. 

In this configuration, unconfined atomizer flow with coflow was directed vertically downward within a 
457 mm^ wire mesh screen. The entire test assembly was surrounded by a flexible plastic enclosure. 

Data were obtained at six axial locations 25, 35, 50, 75, 150, and 300 mm from the exit plane of the injector. 
The details of the injector geometry are shown in Figure 6.6-1. 

The set of governing partial differential equations applied to swirling flow consists of equations for conti- 
nuity, axial, radial, and azimuthal momenta, e, and six Reynolds stress components. The finite volume 
approach (Patankar, 1980) is used to reduce the continuous equations to a set of coupled discrete equa- 
tions. The numerical solutions are obtained using the flux-spline differencing scheme (Varejao, 1979). 

A calculation procedure for elliptic flow requires boundary conditions on all boundaries of the computa- 
tional domain. Four kinds of boundaries need consideration, namely inlet, axis of symmetry, outlet, and 
entrainment. At the inlet boundary, which was located at the first measurement plane, the measured pro- 
files of U, V, W, and the Reynolds stresses were prescribed. The k profile was obtained from the mea- 
sured normal stresses (Figure 6.6-2). The kinetic energy distribution and the measured shear stress pro- 
file were used to derive the inlet dissipation rate values (Equation 69). 

At the axis of symmetry, the radial velocity, shear stresses, and the radial gradients of other variables are 
set to zero. All streamwise gradients are presumed zero in the exit plane of the calculation domain. 

Along the entrainment boundary, which was placed sufficiently far from the axis of symmetry, the quan- 
tity rV was assumed constant. In addition, the axial velocity U was assumed zero, and k and e were as- 
signed arbitraily low values yielding an eddy viscosity, jit = lop. 

The coupled equations and boundary conditions are solved numerically in a sequential manner using the 
staggered grids for velocities and shear stresses. The main advantage of staggering the locations of 
stresses is to enhance the numerical stability, a result of high coupling between the shear stresses and the 
related mean strains. The iteration sequence employs the SIMPLER algorithm (Patankar, 1980) to handle 
the coupling between the continuity and momentum equations. The algebraic equations are solved using 
a line-by-line TDM A. 

A nonuniform grid of 61 x 57 in the axial and radial directions is used for computations (Figure 6.6-3). A 
finer grid spacing is assigned near the inlet, centerline, and in the shear layer. The tabulated axial (x) and 
radial (r) grid points are presented in Table 6.6-1. The computational domain extends from the first exper- 
imental location at x = 25 mm to 450 mm downstream of the jet exit. Since the measured flow does not 
show any x dependence at x > 300 mm, the specified condition = 0 at the exit plane of calculation do- 
main is reasonable. 

The presented calculations have been made using the flux-spline differencing scheme and the calculated 
results are essentially free of false diffusion. Therefore, the discrepancy between the experimental data 
and the prediction can be related to either the improper inlet boundary conditions or the deficiency of the 
turbulence model. As regards the inlet conditions, all quantities except the dissipation rate (e) were pre- 
scribed from the measurement. The uncertainties in the derivation of the inlet e profile would adversely 
affect the calculation at downstream locations. 

Predictions of mean and turbulence fields obtained from the DSM are discussed first. Comparison of the 
calculated mean axial velocity profiles with the experimental data is presented in Figure 6.6-4. In this 
case, the DSM closure mimics the data trend reasonably well. The DSM closure prediction of axial veloc- 
ity is in relatively good agreement with the data. The calculated profiles are in accordance with the data 
trend, although the centerline velocity values are underpredicted. The differences between the model 
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Table 6.6-1. 

Airblast injector flow grid definition. 


STREAMWISE COORDINATES OF THE GRID 


I DX X XU 


1 

0 . 000E+00 

0 . 000E + 00 

0 . OOOE + OO 

2 

7 . 750E-04 

7.750E-04 

0 . 000E+00 

3 

1 . 625E-03 

2 . 400E-03 

1 . 550E-03 

4 

2 . 050E-03 

4 . 450E-03 

3.250E-03 

5 

2.500E-03 

6 . 950E-03 

5 . 650E-03 

6 

2.700E-03 

9 . 650E-03 

8.250E-03 

7 

2.937E-03 

1 . 259E-02 

1 . 105E-02 

8 

3.313E-03 

1 . 590E-02 

1 . 4 1 2E-02 

9 

3.888E-03 

1 . 979E-02 

1 . 768E-02 

10 

4 . 587E-03 

2.437E-02 

2. 190E-02 

11 

5 . 250E-03 

2.962E-02 

2.685E-02 

12 

5.775E-03 

3 . 540E-02 

3 . 240E-02 

13 

6 . 238E-03 

4 . 164E-02 

3.840E-02 

14 

6.71 2E-03 

4 . 835E-02 

4.487E-02 

15 

7 . 062E-03 

5. 54 IE-02 

5 . 1 82E-02 

16 

7.200E-03 

6 . 26 IE-02 

5 . 900E-02 

17 

7.250E-03 

6.986E-02 

6 . 622E-02 

18 

7.288E-03 

7.715E-02 

7.350E-02 

19 

7.300E-03 

8.445E-02 

8.080E-02 

20 

7 . 300E-03 

9. 175E-02 

8.810E-02 

21 

7 . 288E-03 

9 . 904E-02 

9.540E-02 

22 

7.275E-03 

1 . 063E-0 1 

1 .027E-01 

23 

7.288E-03 

1 . 136E-01 

1 .099E-01 

24 

7.300E-03 

1 . 2O9E-01 

1 . 172E-01 

25 

7.300E-03 

1 . 282E-0 1 

1 .245E-01 

26 

7.3G0E-03 

1 . 355E-0 1 

1 .318E-01 

27 

7 . 300E-03 

1 . 428E-0 1 

1 . 391E- 0 1 

28 

7.300E-03 

1 . 50 1E-0 1 

1 . 464E-0 1 

29 

7 . 300E-03 

1 . 574E-01 

1 .537E-01 

30 

7.300E-03 

1 . 647E-01 

1 .610E-01 

31 

7.300E-03 

1 . 720E-01 

1 . 683E-01 

32 

7 . 300E-03 

1 . 793E-0 1 

1 . 756E-0 1 

33 

7 . 300E-03 

1 . 866E-01 

1 .829E-01 

34 

7.300E-03 

1 . 939E-01 

1 .902E-01 

35 

7.300E-03 

2.01 2E-01 

1 .975E-01 

36 

7.300E-03 

2.085E-01 

2 . 048E-0 1 

37 

7.300E-03 

2. 158E-01 

2. 121E-01 

38 

7 .300E-03 

2 . 23 1E-0 1 

2.194E-01 

39 

7.200E-03 

2.303E-01 

2 . 267E-0 1 

40 

7.200E-03 

2.375E-01 

2.338E-01 

41 

7.400E-03 

2.449E-01 

2.41 IE-01 

42 

7.400E-03 

2.523E-01 

2.486E-01 

43 

7.325E-03 

2.596E-01 

2.559E-01 

44 

7.325E-03 

2.669E-01 

2 . 633E-0 1 

45 

7 , 275E-03 

2.742E-01 

2 . 706E-0 1 

46 

7 . 275E-03 

2.815E-01 

2.778E-01 

47 

7.300E-03 

2 . 888E-0 1 

2.851E-01 
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Table 6.6-1. 

Airblast injector flow grid definition (cont). 


48 

7.450E-03 

2.962E-01 

2.925E-01 

49 

7.775E-03 

3.040E-01 

3.001E-01 

50 

7.975E-03 

3.120E-01 

3.080E-01 

51 

8 .500E-03 

3.205E-01 

3. 160E-01 

52 

9.500E-03 

3.300E-01 

3.250E-01 

53 

1 .000E-02 

3.400E-01 

3.350E-01 

54 

1 .000E-02 

3.500E-01 

3.450E-01 

55 

1 . 125E-02 

3.612E-01 

3.550E-01 

56 

1 . 375E-02 

3 . 750E-0 1 

3.675E-01 

57 

1 . 500E-02 

3.900E-01 

3.825E-01 

58 

1 . 500E-02 

4.050E-01 

3.975E-01 

59 

1 .500E-02 

4.200E-01 

4. 125E-01 

60 

1 .500E-02 

4 . 350E-0 1 

4.275E-01 

61 

7.500E-03 

4 .425E-01 

4.425E-01 


TRANSVERS 

COORDINATES OF 

THE GRID 


J 

DY 

Y 

YV 

1 

0 . 000E+00 

5.000E-04 

0 . 000E+00 

2 

2.500E-04 

7.500E-04 

5 . 000E-04 

3 

7 . 500E-04 

1 .500E-03 

1 .000E-03 

4 

1 . OOOE-03 

2 . 500E-03 

2. 000E-03 

5 

1 . 000E-03 

3 .500E-03 

3. 000E-03 

6 

1 .000E-03 

4 .500E-03 

4. 000E-03 

7 

1 .OOOE-03 

5.500E-03 

5. 000E-03 

8 

1 . 000E-03 

6.500E-03 

6. 000E-03 

9 

1 . 000E-03 

7.500E-03 

7. 000E-03 

10 

1 . 000E-03 

8 .500E-03 

8. OOOE-03 

1 1 

1 . 000E-03 

9.500E-03 

9. 000E-03 

12 

1 .000E-03 

1 .050E-02 

1 . 000E-02 

13 

1 .000E-03 

1 . 150E-02 

1 . 100E-02 

14 

1 . 000E-03 

1 .250E-02 

1 .200E-02 

15 

1 . 000E-03 

1 . 350E-02 

1 .300E-02 

16 

1 .000E-O3 

1 . 450E-02 

1 . 400E-02 

17 

1 .000E-03 

1 . 550E-02 

1 . 500E-02 

18 

1 .000E-03 

1 .650E-02 

1 .600E-02 

19 

1 . 000E-03 

1 . 750E-02 

1 .700E-02 

20 

1 .O00E-03 

1 .850E-02 

1 . 800E-02 

21 

1 . 000E-03 

1 .950E-02 

1 . 900E-02 

22 

1 . O00E-03 

2.050E-02 

2.000E-02 

23 

1 . 000E-03 

2. 150E-02 

2. 100E-02 

24 

1 . 000E-03 

2.250E-02 

2.200E-02 

25 

1 .000E-03 

2.350E-02 

2.300E-02 

26 

1 .000E-03 

2.450E-02 

2.400E-02 

27 

1 . 000E-03 

2.550E-02 

2.500E-02 

28 

1 . 000E-03 

2.650E-02 

2 . 600E-02 

29 

1 .000E-03 

2.750E-02 

2.700E-02 

30 

1 .OOOE-03 

2.850E-02 

2.800E-02 

31 

1 .000E-03 

2.950E-02 

2.900E-02 
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Table 6.6-1. 

Airblast injector flow grid definition (cont). 


32 

1 . 000E-03 

3.050E-02 

3.000E-02 

33 

1 .000E-03 

3. 150E-02 

3. 100E-02 

34 

1 .OOOE-03 

3.250E-02 

3.200E-02 

35 

1 .500E-03 

3.400E-02 

3.300E-02 

36 

2. OOOE-03 

3 . 600E-02 

3.500E-02 

37 

2. 000E-03 

3.800E-02 

3.700E-02 

38 

2. 000E-03 

4.000E-02 

3.900E-02 

39 

2. 000E-03 

4.200E-02 

4. 100E-02 

40 

2. 000E-03 

4 .400E-02 

4.300E-02 

41 

2. 000E-03 

4 .600E-02 

4.500E-02 

42 

2. 000E-03 

4.800E-02 

4.700E-02 

43 

2. OOOE-03 

5.000E-02 

4.900E-02 

44 

2. OOOE-03 

5.200E-02 

5. 100E-02 

45 

2. 000E-03 

5.400E-02 

5.300E-02 

46 

2. OOOE-03 

5.600E-02 

5.500E-02 

47 

2. 000E-03 

5.800E-02 

5.700E-02 

48 

2. 000E-03 

6.000E-02 

5.90DE-02 

49 

2 . 000E-03 

6.200E-02 

6. 100E-02 

50 

2. 000E-03 

6 . 400E-02 

6.300E-02 

51 

2. 000E-03 

6.600E-02 

6.500E-02 

52 

2. 000E-03 

6 .800E-02 

6.700E-02 

53 

2. 000E-03 

7.000E-02 

6.900E-02 

54 

2. OOOE-03 

7.200E-02 

7. 100E-02 

55 

2. 000E-03 

7 . 400E-02 

7.300E-02 

56 

1 .500E-03 

7.550E-02 

7.500E-02 

57 

5.000E-04 

7.600E-02 

7.600E-02 


and data become more severe in the developing region. The discrepancies are due to the insufficient dif- 
fusion process caused by underprediction of radial normal stress component. At downstream, the calcu- 
lated velocity field is in good agreement with data, however some minor discrepancies appear in that 
zone. 


Comparisons of the predicted turbulence intensities from the DSM closure are shown in Figures 6.6-5, 6.6- 
6, and 6.6-7. The differences between the calculated and experimental results are significant in the case of 
the normal stresses. The maximum axial turbulence intensity has been underpredicted in most of the 
flow region, however, the trend was accurately predicted. The same trends are observed for all normal 
stress components. The discrepancies between the model and data clearly indicate the lack of perfor- 
mance of the pressure-strain model. One reason for this could be the lack of an appropriate model for the 
rapid part, which was not considered here. The predicted turbulent shear stress from the DSM has been 
compared with the experimental data in Figure 6.6-8. The positive peak in the shear stress profile corre- 
sponds to the shear layer between the two streams. The agreement between the calculation and the ex- 
perimental data is relatively good, although the peak values are not well predicted. 

The k-e turbulence model was also applied for the same configuration to establish the suitability of the 
model. The predicted mean axial velocity, turbulent kinetic energy, and shear profiles at the various axial 
stations are shown in Figures 6.6-9, 6.6-10, and 6.6-11. The predictions of axial velocity profiles are in fair 
agreement with the data. The velocity profiles are relatively flat near the centerline. The maximum veloc- 
ity occurs near the centerline which constantly decreases. The deceleration in axial velocity is due to the 
pressure effects. The agreement between the predicted and experimental values of kinetic energy is as 
good as that for the axial velocity. Even though the trends are similar, the predicted kinetic energy levels 
are smaller than those derived from the measurements. As regards the turbulent shear stress, the calcu- 
lated profiles only qualitatively agree with the experimental data. The predictions are not as good as 
those calculated by the DSM closure. 
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Figure 6.6-2. Airblast injector flow inlet conditions. 
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Figure 6.6-4. Comparison of calculated mean axial velocity profiles by DSM with data. 
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Figure 6.6-5. Comparison of calculated rms axial velocity profiles by DSM with data. 
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Figure 6.6-6. Comparison of calculated rms radial velocity profiles by DSM with data. 




Figure 6.6-7. Comparison of calculated rms tangential velocity profiles by DSM with data. 
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Figure 6.6-8. Comparison of calculated shear stress profiles by DSM with data. 
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Figure 6.6-9. Comparison of calculated mean axial velocity profiles by k-e with data. 
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Figure 6.6-10. Comparison of calculated turbulent kinetic energy profiles by k-e with data. 
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Figure 6.6-11. Comparison of calculated turbulent shear stress profiles by k-e with data. 
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VII. CONCLUDING REMARKS AND RECOMMENDATIONS 


7.1 SUMMARY OF THE PRESENT WORK 

A combined experimental/analytical investigation was conducted to validate conventional and improved 
turbulent and spray models. A systematic experimental program was conducted to collect benchmark 
quality data for the mean and fluctuating quantities in the developing region of single-phase and particle- 
laden flows with two mass loading ratios. Detailed velocity measurements were made using phase/Dop- 
pler equipment. A detailed specification of the flow parameters in the upstream region was provided. 
These are used as inlet conditions to start the computations. The test problems used ranges from simple 
flows to more complex flows to encompass the range of complexities involved in combustion flows. 

These are 

• single round jet 

• single annular jet 

• single swirling annular jet 

• coaxial jets 

• coaxial jets with swirling anr ular flow 

• airblast injector 

Closure of the Reynolds equation was achieved by three different levels of models: k-e, algebraic second- 
moment (ASM), and differential second-moment (DSM) closure . Effects of pressure-strain correlation 
models on these closure were also investigated. For turbulent particle-laden jet flow, the closure was af- 
fected by extending constant-density single-phase turbulence models to the continuous (carrier) phase of 
two-phase flow. The dispersed (discrete) phase was treated by both deterministic and stochastic-Lagran- 
gian technique. 

To reduce the effect of numerical (false) diffusion on the predicted results, the linear flux-spline scheme 
was used to solve several two-dimensional flows. For a given number of grid points, the flux-spline 
scheme produces results that are superior to those from tine (lower-order) power-law differencing 
scheme. In addition, it has the potential of providing a grid independent solution without requiring an 
excessive number of grid points. 

Models employing transport equations for the individual turbulent stress components simulate the 
turbulent processes more realistically and are, therefore, potentially more general compared to the 
simpler models. It should be noted that the effectiveness of turbulence model predictions could be 
obscured to some extent by competing factors: boundary conditions, oscillatory phenomena, and 
numerical diffusion. A significant contribution from any of the aforementioned factors tends to 
invalidate conclusions regarding the superiority or inferiority of a given turbulence model. The 
turbulence model cannot compensate for inadequacy in this area. 

7.2 RECOMMENDATION 

Although some of the models described in this study, and in particular the k-e model, have been shown 
to work well in many situations, there is much room for further development. The e-equation in its pre- 
sent form appears not to be sufficiently universal and should be improved. As observed by many, this 
equation is responsible for performance of most models. Ideas to use several length-scale equations for 
different directions or different processes are promising and should be developed further. 

The model assumptions for the pressure-strain correlation are also not very satisfactory and need im- 
provement. Proposals for the behavior of <&iji in homogeneous flows have only gone further than that of 
Rotta's linear model by including further terms in a series expansion about the isotropic homogeneous 
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state. However, optimization of the coefficients of the terms in the expansion on the basis of available ex- 
perimental data is a very difficult task indeed. 

It seems unlikely that any serious proposal for will be made in the near future, and emphasis should 
be placed on developing a better approximation for 

The derivation and validation of higher-order closure schemes holds the greatest potential for turbulence 
model improvement for strongly swirling flows. Efforts to find a stable and higher order (order of terms 
retained in an equivalent Taylor series expansion) differencing scheme thai can eliminate numerical dif- 
fusion should continue. This is especially important in the case of two-phase and/ or reacting flows be- 
cause of the coupled nonlinearities which exist between the chemical and fluid mechanical processes. 

An intensive submodel validation and development efforts, especially for the ASM closure, DSM closure, 
and stochastic spray model in conjunction with probability distribution function (PDF) approach should 
be continued. 

Unfortunately, in many instances there is a lack of quality data relevant to gas turbine combustion. Many 
modeling assumptions are similar to the constant-density, Reynolds t tress closure. Therefore, further ex- 
periments with more emphasis on turbulent scalar fluxes and density correlations are needed to support 
or to improve these assumptions. 
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TURBULENT FLOW EQUATIONS FOR THE k-e MODEL 

The transport equations presented in the previous sections reduce in axisymmetric coordinates (x, r) to 
the following: 

Continuity Equation 
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Turbulence Model 


Ht = C u P 
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TURBULENT FLOW EQUATIONS FOR DSM CLOSURE 


The modeled Reynolds stress transport equations presented in the tensor notation reduce in axisymmetric 
coordinates (x, r) to the following equations. 
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Equation for w : 
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